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Abstract. These lecture notes review the foundations and some applications of light- 
cone quantization. First I explain how to choose a time in special relativity. Inclusion 
of Poincare invariance naturally leads to Dirac's forms of relativistic dynamics. Among 
these, the front form, being the basis for light -cone quantization, is my main focus. 
I explain a few of its peculiar features such as boost and Galilei invariance or sep- 
aration of relative and center-of-mass motion. Combining light-cone dynamics and 
field quantization results in light -cone quantum field theory. As the latter represents 
a first-order system, quantization is somewhat nonstandard. I address this issue us- 
ing Schwinger's quantum action principle, the method of Faddeev and Jackiw, and 
the functional Schrodinger picture. A finite-volume formulation, discretized light -cone 
quantization, is analysed in detail. I point out some problems with causality, which are 
absent in infinite volume. Finally, the triviality of the light -cone vacuum is established. 
Coming to applications, I introduce the notion of light -cone wave functions as the so- 
lutions of the light-cone Schrodinger equation. I discuss some examples, among them 
nonrelativistic Coulomb systems and model field theories in two dimensions. Vacuum 
properties (like chiral condensates) are reconstructed from the particle spectrum ob- 
tained by solving the light-cone Schrodinger equation. In a last application, I make 
contact with phenomenology by calculating the pion wave function within the Nambu 
and Jona-Lasinio model. I am thus able to predict a number of observables like the pion 
charge and core radius, the r.m.s. transverse momentum, the pion structure function 
and the pion distribution amplitude. The latter turns out to be the asymptotic one. 



1 Introduction 

The nature of elementary particles calls for a synthesis of relativity and quantum 
mechanics. The necessity of a quantum treatment is quite evident in view of the 
microscopic scales involved which are several orders of magnitude smaller than in 
atomic physics. These very scales, however, also require a relativistic formulation. 
A typical hadronic scale of 1 fm, for instance, corresponds to momenta of the 
order of p ~ hc/lim ~ 200 MeV. For particles with masses M ~ 1 GeV, this 
implies sizable velocities v ~ p/M ~ 0.2 c. 

It turns out that the task of unifying the principles of quantum mechan- 
ics and relativity is not a straightforward one. One can neither simply extend 
ordinary quantum mechanics to include relativistic physics nor quantize rela- 
tivistic mechanics using the ordinary correspondence rules. Nevertheless, Dirac 
and others have succeeded in formulating what is called "relativistic quantum 
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mechanics" , which has become a subject of text books since — see e.g. (Bjorkcn 
and Drell 1964; Yndurain 1996). It should, however, be pointed out that this 
formulation, which is based on the concept of single-particle wave-functions and 
equations, is not really consistent. It does not correctly account for relativistic 
causality (retardation effects etc.) and the existence of antiparticles. As a result, 
one has to struggle with issues like the Klein paradox 1 , the definition of position 
operators (Newton and Wigner 1949) and the like. 

The well-known solution to these problems is provided by quantum field the- 
ory, with an inherently correct description of antiparticles that entails relativistic 
causality. In contrast to single-particle wave mechanics, quantum field theory is 
a (relativistic) many body formulation that necessarily involves (anti-)particle 
creation and annihilation. Physical particle states are typically a superposition 
of an infinite number of 'bare' states, as any particle has a finite probability 
to emit or absorb other particles at any moment of time. A pion, for example, 
would be represented in terms of the following Fock expansion, 



where the tj) n are the probability amplitudes to find n particles (quarks q, anti- 
quarks q or gluons g) in the pion. With the advent of QCD, however, a concep- 
tual difficulty concerning this many-particle picture has appeared. At low energy 
or momentum transfer, hadrons, the bound states of QCD, are reasonably de- 
scribed in terms of two or three constituent quarks and thus as few-body systems. 
These 'effective' quarks Q are dressed so that they gain an effective mass of the 
order of 300-400 MeV. They are used as the basic degrees of freedom in the 
'constituent quark model'. This model yields a reasonable mass spectroscopy 
of hadrons (Lucha et al. 1991; Mukherjee et al. 1993), but its foundations are 
not very well established theoretically. First, a nonrelativistic treatment of light 
hadrons is not justified (see above). Second, the model violates many symme- 
tries of QCD (in particular chiral symmetry). Third, it is rather unclear how a 
constituent picture can arise in a quantum field theory such as QCD. 

In principle, in order to confirm the constituent quark model, one would have 
to solve the 'QCD Schrodinger equation' for hadron states |hadron) of mass Mh, 



and check whether the eigenstates are reasonably well described in terms of the 
constituent valence states \QQ) or \QQQ). This is a very hard problem. A more 
moderate goal would be to 'relativize' the constituent quark model, ideally in 
such a way that it respects the symmetries of QCD. I will discuss this attempt 
in detail at the end of these lectures. 

To arrive at this point, there is, of course, some way to go. Let me start with 
the following claim. A particularly useful approach for our purposes is based 
on a somewhat unorthodox choice of the 'time arrow' within special relativity: 

1 For a nice recent discussion, see (Holstein 1998). 



k) = ip2\qq) + ip3\qqg) + ^\qqqq) + ■■■ 



(l) 



Hqcd I hadron) = Mh| hadron) , 



(2) 
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instead of the ordinary 'Galileian time' t, I choose 'light-cone time', x + = t — z/c. 
In the course of these lectures 2 , this claim will be substantiated step by step. 

I will begin with some general remarks on relativistic dynamics (Section 2). 
As a paradigm example I discuss the free relativistic particle which is the proto- 
type of a reparametrization invariant system. I show that the choice of the time 
parameter is not unique as it corresponds to a gauge fixing, the purpose of which 
is to get rid of the reparametrization redundancies. By considering the stabil- 
ity subgroups of the Poincare group, one finds that there are essentially three 
reasonable choices of 'time' for a relativistic system, corresponding to Dirac's 
'instant', 'point' and 'front' form, respectively. The latter choice is the basis of 
light-cone dynamics, the main features of which will be discussed in the last 
part of Section 2. 

Section 3 is devoted to light-cone field quantization. I show how the Poincare 
generators are defined in this case and utilize Schwinger's quantum action prin- 
ciple to derive the canonical commutators. This is the first method of quanti- 
zation to be discussed. The relation between equal-time commutators, the field 
equations and their solutions for different initial and/or boundary conditions is 
clarified. 

It turns out that light-cone field theories, being of first order in the velocity, 
generally are constrained systems which require a special treatment. I rederive 
the canonical light-cone commutators using a second method of quantization 
(based on phase space reduction) due to Faddeev and Jackiw. I extend this 
discussion to light cone quantization in finite volume and point out possible 
problems with causality in this approach. Going back to infinite volume, I in- 
troduce a third method of quantization, the functional Schrodingcr picture, and 
combine it with the light-cone formalism. I close this section with a discussion 
of the presumably most spectacular feature of light-cone quantum field theory, 
the triviality of the vacuum. 

As a prelude to the applications I introduce the notion of light-cone wave 
functions in Section 4. I show how light-cone wave functions can be obtained 
by solving the light-cone Schrodingcr equation. As examples, I discuss nonrela- 
tivistic wave functions as they occur in hydrogen-like systems, some model field 
theory in 1+1 dimensions and a simple Gaussian model. 

In Section 5, 1 finally make contact with phenomenology. I calculate the light- 
cone wave function of the pion within the Nambu and Jona-Lasinio model. This 
model is known to provide a good description of spontaneous chiral symmetry 
breaking, as it is governed by the same symmetry group as low-energy QCD. 
With the pion wave function at hand I derive a number of observables like the 
pion charge and core radius, the electromagnetic form factor, the r.m.s. transverse 
momentum and the pion structure function. I conclude with a calculation of the 
pion distribution amplitude. 



2 Naturally, these lectures cannot cover all aspects of light-cone field theory. For a 
comprehensive recent review on the subject see (Brodsky, Pauli, and Pinsky 1998). 



4 Thomas Heinzl 



2 Relativistic Particle Dynamics 

The physically motivated desire to describe hadrons as bound states of a small, 
fixed number of constituents is our rationale to go back and reanalyze the relation 
between Hamiltonian quantum mechanics and relativistic quantum field theory. 

Quite generally, bound states are obtained by solving the Schrodinger equa- 
tion, 

in-j^Mr)) = hw{t)) , (3) 

for normalized, stationary states, 

|V>(T))=e-^|V(0)>. (4) 
This leads to the bound-state equation 

mm) = E\m) , (5) 

where E is the bound state energy. We would like to make this Hamiltonian 
formalism consistent with the requirements of relativity. It is, however, obvious 
from the outset that this procedure will not be manifestly covariant as it singles 
out a time r (and an energy E, respectively). Furthermore, it is not even clear 
what the time r really is as it does not have an invariant meaning. 



2.1 The Free Relativistic Point Particle 

To see what is involved it is sufficient to consider the classical dynamics of a 
free relativistic particle. We want to find the associated canonical formulation 
as a basis for subsequent quantization. We will proceed by analogy with the 
treatment of classical free strings which is described in a number of textbooks 
(Green et al. 1987; Polyakov 1987). Accordingly, the relativistic point particle 
may be viewed as an infinitely short string. 

Recall that the action for a relativistic particle is essentially given by the arc 
length of its trajectory 

S = —ms\2 = — m J ds . (6) 

This action 3 is a Lorentz scalar as 

ds = y/g^xW (7) 
is the (infinitesimal) invariant distance. We can rewrite the action (6) as 

S=-mJ^ ds y/x^ = J dsL(s), (8) 
3 We work in natural units, % = c = 1. 
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in order to introduce a Lagrangian L(s) and the four velocity x M = dx^/ds. The 
latter obeys 

x 2 = x^ = 1 , (9) 



as the arc length provides a natural parametrization. Thus, i M is a time-like 
vector, and we assume in addition that it points into the future, x° > 0. In 
this way we guarantee relativistic causality ensuring that a real particle passing 
through a point P will always propagate into the future light cone based at P. 

We proceed with the canonical formalism by calculating the canonical mo- 
menta as 

f)T 

pf* = -— = mi" . (10) 
oxf, 

These are not independent, as can be seen by calculating the square using (9), 

p 2 = m 2 x 2 = m 2 , (11) 

which, of course, is the usual mass-shell constraint. This constraint indicates 
that the Lagrangian L(s) defined in (8) is singular, so that its Hessian W^ u with 
respect to the velocities, 

^ - 7^4- = I iT -^-)= -mGT - srr) . (12) 



is degenerate. It has a zero mode given by the velocity itself, 

W^x v = . (13) 

The Lagrangian being singular implies that the velocities cannot be uniquely 
expressed in terms of the canonical momenta. This, however, is not obvious 
from (10), as we can easily solve for the velocities, 

x»= p v/m. (14) 

But if one now calculates the canonical Hamiltonian, 

H c = -p^ — L = -mx 2 + mx 2 = , (15) 

one finds that it is vanishing! It therefore seems that we do not have a generator 
for the time evolution of our dynamical system. In the following, we will analyze 
the reasons for this peculiar finding. 

First of all we note that the Lagrangian is homogeneous of first degree in the 
velocity, 

L(ax") = aL{x^) . (16) 
Thus, under a reparamctrization of the world-line, 



shW, x"(s) ^x»(s'(s)) , 



(17) 
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where the mapping s s' is one-to-one with ds' /ds > (orientation conserv- 
ing), the Lagrangian changes according to 

L(dx"/ds) = L^(dx"/ds')(ds'/ds)^ = (ds'/ds)L(dx> 1 /ds') . (18) 

This is sufficient to guarantee that the action is invariant under (17), that is, 
reparametrization invariant, 



f S2 r s ' 2 ds ds' 

/ dsLidx^/ds) = / ds' -L^—L(dx^/ds') = S' 
J Sl J s , ds' ds 



(19) 



if the endpoints remain unchanged, s\ : 2 = s' 12 - On the other hand, L is homo- 
geneous of the first degree if and only if Euler's formula holds, namely 

L = dx~^ = ■ (20) 

This is exactly the statement (15), the vanishing of the Hamiltonian. Further- 
more, if wc differentiate (20) with respect to x^, we recover (13) expressing 
the singular nature of the Lagrangian. Summarizing, we have found the general 
result (Hanson et al. 1976; Sundermeyer 1982) that if a Lagrangian is homoge- 
neous of degree one in the velocities, the action is reparametrization invariant, 
and the Hamiltonian vanishes. In this case, the momenta are homogeneous of 
degree zero, which renders the Lagrangian singular. 

The reparametrization invariance is generated by the first class constraint 
(Dirac 1964; Sundermeyer 1982), 

9 = p 2 -to 2 = , (21) 

as can be seen as follows. From the canonical one form —g^p^dx" we read off 
the Poisson bracket 

{x» , p v ] = -(T , (22) 
and calculate the change of the coordinate 

8x^ = {x» , 95e} = -2p^8e = -2mx»5e = x^St 

= x^(t + St) - x»(t) = .x^(r') - x»(t) . (23) 

Thus, the reparametrization (17) is indeed generated by the constraint (21). 

Reparametrization invariance can be viewed as a gauge or redundancy sym- 
metry. The redundancy consists in the fact that a single trajectory (world-line) 
can be described by an infinite number of different parametrizations. The physi- 
cal objects, the trajectories, are therefore equivalence classes obtained by identi- 
fying ('dividing out') all reparametrizations. The method to do so is well known, 
namely gauge fixing. For the case at hand, this corresponds to a particular 
choice of parametrization, or, more physically, to the choice of a time parame- 
ter t. This amounts to choosing a foliation of space-time into space and time. 
Minkowski space is thus decomposed into hypersurfaces of equal time, t = const, 
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which in general are three-dimensional objects, and the time direction 'orthogo- 
nal' to them. The time development thus continuously evolves the hypersurface 
U : t = t into S\ : t = t\ > t . Put differently, initial conditions provided 
on E together with the dynamical equations (being differential equations in r) 
determine the state of the dynamical system on Si . 

Practically, the (3+l)-foliation is done as follows. We introduce some arbi- 
trary coordinates, = £ Q (a;), which may be curvilinear. We imagine that three 
of these, say £ l , i = 1,2,3, parametrize the three-dimensional hypersurface S, 
so that the remaining one, £°, represents the time variable, i.e. r = £,°(x). This 
equation can equivalently be viewed as a gauge fixing condition. 

The first question to be addressed is: what is a 'good' choice of time? There 
are two criteria to be met, namely existence and uniqueness. Existence means 
that the equal time hypersurface S should intersect any possible world-line, 
while uniqueness requires that it docs so once and only once. Mathematically, 
uniqueness can be analysed in terms of the Faddeev-Popov (FP) 'operator', 
which is given by the Poisson bracket of the gauge fixing condition with the 
constraint (evaluated on £), 

FP = {?{x) , 9} = {e(x) , p 2 } = -2 = -2 N ■ p . (24) 

Here, we have introduced the normal N on E, 

d£°{x) 



iV"(ar) 



dx^ 



(25) 



which will be important later on. The statement now is that uniqueness is 
achieved (for a single degree of freedom) if the FP operator does not vanish, 
i.e. if N ■ p ^ 0. Generically, this means that the particle trajectory must not be 
parallel to the hypersurface S of equal time. 

As an aside we remark that this is completely analogous to the reasoning in 
standard gauge (field) theory. There, the constraint 8 is given by Gauss's law 
which generates gauge transformations A — > A + Duo, D denoting the covariant 
derivative. For a gauge fixing x[A] = 0, the equation corresponding to (24) 
becomes 

where all (functional) derivatives are to be evaluated on the gauge fixing hyper- 
surface, x = 0. 

Let us now perform an analysis of the canonical formalism for a general choice 
of hypersurface S. For this we need some notation. We write the line element as 

ds 2 = g^dx^dx" - g ^ — —dCde = h afj {0d(, a de . (27) 

Introducing a vierbein e^(£), the metric h a f}{^) is alternatively given by 
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The transformation x —* £ is well known from general relativity, where it cor- 
responds to the transformation from a local inertial frame described by the flat 
metric to a nonincrtial frame with coordinate dependent metric h a p(t;). For 
our purposes we write this metric in a (3+l)-notation as follows, 

h ^-[h i0 h ij ) = [h -h) ■ (29) 

Of particular interest is the component hoo, which explicitly reads 

hoo = = fl^Vo = n , (30) 

where we have defined the unit vector in ^"-direction 

n» = ^ = e%=i» % (31) 
which thus is the new four-velocity. It is related to the normal vector via 

u.JV^eV; = I^S£ = 1- (32) 
The normal vector N enters the inverse metric which we write as follows, 

\9io 9ij ) V 9 -G ) 

The hij are the metric components associated with the hypersurface. The in- 
variant distance element (27) thus becomes (with h oi = hi), 

ds 2 = hoodt d£ + 2h 0i dt°dC + hijdgd? , 

n 2 + 2h i ^f+h ij ^f^f)dT 2 =h(r)dT 2 , (34) 
dr dr dr J 

where, in the second step, we have used that £° = r. In the last identity we have 
defined a world-line metric or cinbcin 

h(r) = x 2 = h a0 i a ^ , (35) 

which expresses the arbitrariness in choosing a time by providing an (arbitrary) 
'scale' for the velocity. Introducing the velocities expressed in the new coordi- 
nates, w l = d^/dr, the world-line metric can be written as 

h{r) = n 2 + 2h,w l + h i3 w i w j . (36) 

Let us develop a canonical formalism for a general choice of the einbein h(r) 
corresponding to the gauge fixing r = £°(x) = 0. The Lagrangian becomes 

L(t) = -m^Uj) , (37) 
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leading to the canonical momenta 

ir a = -^ = ^h af3 i a =e£p„. (38) 

We see that the einbein h is appearing all over the place. The canonical Hamil- 
tonian is expressed in terms of the inverse metric h a ^ , 

H c = -n a i a -L = -^-(h^n^p - m 2 ) = . (39) 

It vanishes (as it should) as it is proportional to the constraint, 

6 = h aP K a K fi -m 2 = p 2 -to 2 =0 . (40) 

The FP operator also depends on the entries of the inverse metric (33), in par- 
ticular the normal vector N, 

FP = N 2 ir + sVj . (41) 

After gauge fixing, the generator of r-evolution, H T = 7r , is obtained by solving 
the constraint (40) for n which assumes the explicit form, 

N 2 ir 2 + 2 3 V 4 7T - G^TTiTTj — to 2 = . (42) 

Depending on the value of A^ 2 , we thus have to consider two different cases. 
The generic one is that the normal A~ on S is time-like, N 2 > 0. In this case, 
the mass-shell constraint is of second order in ttq, so that there are two distinct 
solutions, 

= ^ {-(g, tt) ± y/(g,7r)*+N*[(Tr,Gir)+m?}} . (43) 

Not unexpectedly, the 'problem' of two different signs in front of the square root 
arises (Gitman and Tyutin 1990). Within quantum mechanics, this is somewhat 
difficult to interpret. Upon 'second quantization', i.e. in the context of quantum 
field theory, one has, of course, the natural explanation in terms of antiparticles. 
As we will not quantize the relativistic point particle, the sign 'problem' is of 
no concern to us. A possible arbitrariness will be removed ad hoc by demanding 
7r > 0. With this additional condition the FP operator becomes 

FP = -2v/( S ,7r) 2 + N 2 [(tt,Gtt) +m 2 } , (44) 

which is clearly nonvanishing for a massive particle, m ^ 0. A gauge fixing with 
A^ 2 > is thus unique and leads to a well-defined description of the T-evolution. 

The second case to be considered is in a sense degenerate. It corresponds to 
a light-like normal, A^ 2 = 0. In this case, the constraint (42) is only of first order 
in -no leading to a single solution, 

(tt , Gtt) + m 2 

n = ± 7 -L^ . 45 

9,tt 
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As a result, there is no 'sign problem' and no 'ugly' square root. Conservation 
of difficulties, however, is at work, because it is no longer obvious whether the 
FP operator, 

FP = -2(<7,7r) , (46) 

is different from zero. Clearly, this is absolutely necessary for (45) to represent 
a well-defined solution. 

At this point, it should be mentioned that the results (43) and (45) are not 
yet the full story. The entries of the inverse metric, A 2 , g and G should actually 
be expressed in terms of the quantities n 2 , h and H defining the induced metric 
on S. So far, it is also not completely clear which choices of these parameters 
actually make sense physically. Of course, the normal A should not be space- 
like as this would imply that 27 contains time-like directions and thus possible 
particle trajectories. In the next subsection I will give some criteria for reasonable 
choices of time. 

Before we come to that let us apply the general formalism to the standard 
choice of 'Galileian' time, r = £°(x) = x° = t. In this case, the surface S : t = 
is an entirely space-like hyperplanc with constant normal vector A — (1, 0) = n. 
The other metric entries are h = g = and H = G = 1. The world-line metric 
(35) thus becomes 

h{t) = x 2 = l-v 2 = l/ 1 2 , (47) 

where 7 is the usual Lorentz contraction factor. The Hamiltonian is obtained in 
line with the second-order case above, 

H t =N-p = p° = Vp 2 + m 2 ~ FP . (48) 

It generates the dynamics via the basic Poisson bracket {x l , p- 7 } = <5 y leading 
to 

x i = {x i ,H t }=p i /p° , (49) 

with p° given by (48). Note that a well-defined time evolution requires a nonva- 
nishing FP operator (which is proportional to p°). 

As already announced, we will discuss alternatives to this standard choice of 
time in the next subsection. 



2.2 Dirac's Forms of Relativistic Dynamics 

To address this issue it is not sufficient to consider only the T-development and 
the associated generator of time translations (i.e. the Hamiltonian) H T . Instead, 
one has to refer to the full Poincarc group to be able to guarantee full relativistic 
invariance. The generators of the Poincare group are the four momenta P M and 
the six operators M^ v which combine angular momenta and boosts according 
to 

V = ^e ijk M jk , (50) 
K l = M 0i , (51) 
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with i, j, k = 1,2,3. These generators are elements of the Poincare algebra which 
is defined by the Poisson bracket relations, 

{P p ,P V } = 0, 

{M" v , P p } = g vp P p - g pp P v , (52) 
, M pa } = g pa M vp ~ g^M 1 "* - g ua M pp + g vp M iia . 

It is well known that the momenta P M generate space-time translations and 
the M pu rotations and Lorentz boosts, cf. (50, 51). In the following we will 
only consider proper and orthochronous Lorentz transformations, i.e. we exclude 
space and time reflections. 

Any Poincare invariant dynamical theory describing e.g. the interaction of 
particles should provide a particular realization of the Poincare algebra. For 
this purpose, the Poincare generators are constructed out of the fundamental 
dynamical variables like positions, momenta, spins etc. An elementary realization 
of (52) is given as follows. Choose the space-time point x^ and its conjugate 
momentum p p as canonical variables, i.e. adopt (22), 

{x» , p"} - -<T ■ (53) 

The Poincare generators are then found to be 

P" =p f * , M pu = x»p v - x v p» , (54) 

as is easily confirmed by checking (52) using (53). An infinitesimal Poincare 
transformation is thus generated by 

SG = -\8u) lxv Mi xv + 6a lx P> x (55) 

in the following way, 

Sx" = {x" , 5G} = Su^Xu + Sa" , 5uj^ = -5uj up . (56) 
The action of the Poincare group on some scalar function F(x) is thus given by 

SF = {F, SG} = d p F 6a„ - \ {x p d v - x v d»)F <Jw M „ . (57) 

Though the realization (54) is covariant, it has several shortcomings. It docs not 
describe any interaction; for several particles the generators are simply the sum 
of the single particle generators. This point, however, is of minor importance to 
us, and will only be touched upon at the end of this subsection. The solution 
of the problem, as already mentioned in the introduction, is the framework of 
local quantum field theory. More importantly, the representation (54) does not 
take into account the mass-shell constraint, p 2 = m 2 , which we already know 
to guarantee relativistic causality as it generates the dynamics upon solving for 
H T . 

To remedy the situation we proceed as before by choosing a time variable 
r, i.e. a foliation of space-time into essentially space-like hypersurfaces £ with 
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time-like or light-like normals N. We have seen that E should be chosen in 
such a way that it intersects all possible world-lines once and only once (exis- 
tence and uniqueness). Apart from this necessary consistency with causality the 
foliation appears quite arbitrary. However, given a particular foliation one can 
ask the question which of the Poincare generators will leave the hypersurface 
E invariant. The set of all such generators defines a subgroup of the Poincare 
group called the stability group Gs of E. The associated generators are called 
kinematical, the others dynamical. The latter map E onto another hypersurface 
E' and thus involve the development in r. One thus expects that the dynami- 
cal generators will depend on the Hamiltonian (and, therefore, the interaction) 
which, by definition, is a dynamical quantity. 

It is clear, however, that the stability group corresponding to a particular 
foliation will be empty if the associated hypersurface looks very irregular and 
thus does not have a high degree of symmetry. One therefore demands in addition 
that the stability group acts transitively on E: any two points on E can be 
connected by a transformation from Ge- With this additional requirement there 
are exactly five inequivalent classes of hypersurfaces (Leutwyler and Stern 1978) 
which are listed in Table 1. 



Table 1. All possible choices of hypersurfaces E: r — const with transitive action of the 
stability group Ge- d denotes the dimension ofGs, that is, the number of kinematical 
Poincare generators; x± = {x 1 ,x 2 ). 



name 


E 










r 






d 


instant 


x° 


= 








t 






6 


light front 


x° 


+ x 3 


= 






t + x 


3 /c 




7 


hyperboloid 


xl 




= a 2 


> 0, x° 


> 


(t 2 - 


2 / 2 

X /c — 


a 2 /c 2 ) 1 ' 2 


6 


hyperboloid 


2 
X 


2 

- Xl 


2 

= a 


> 0, x° 


> 




2 / 2 


-a 2 /c 2 ) 1/2 


4 


hyperboloid 


Xq 


- x\ 


= a 


> 0, x° 


> 


(t 2 - 


xi/c - 


a 2 /c 2 r 2 


4 



The first three choices have already been found by (Dirac 1949) in his seminal 
paper on 'forms of rclativistic dynamics'. He called the associated forms the 
'instant', 'front' and 'point' forms, respectively. These are the most important 
choices as the remaining two forms have a rather small stability group and thus 
are not very useful. We have only listed them for the sake of completeness. 

It is important to note that for all forms one has limc_>oo r = t, which means 
that in the nonrelativistic case there is only one possible foliation leading to the 
absolute Galileian time t. This is consistent with the fact that there is no limiting 
velocity in this case implying that particle trajectories can have arbitrary slope 
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(tangent vector). Therefore, the hypersurface S nr : t = const is the only one 
intersecting all possible world-lines. For other choices, the criterion of existence 
introduced in the last subsection would be violated. 

To decide which of the Poincare generators are kinematical, we use the general 
formula (57) describing their action. Imagine that £ is given in the form S : 
t = £°(x) = F(x) as in Table 1. If P^ or are kinematical for some u or v, 
then, for these particular superscripts, the components of the gradient and rotor 
of F have to vanish on E, 

d»F = 0, {x^d v - x'd^F = . (58) 

In terms of the normal vector ./V these equations become 

= , aPN" - x u N" = , (59) 

which again will hold for some of the superscripts u and/or v, if E has nontrivial 
stabilizer. The distinction between kinematical and dynamical is thus completely 
encoded in the normal vector N . 

The choice of Galileian time r = t is of course the most common one also in 
the relativistic case, and we have discussed it briefly at the end of the last sub- 
section. To complete this discussion, we construct the associated representation 
of the Poincare generators on E : t = 0. The idea is again to explicitly saturate 
the constraint p 2 = m 2 by solving for H t = p° = N ■ p = (p 2 + m 2 ) 1 / 2 , and 
setting x° = in (54). 

As a result, we obtain the following (3+l)-representation of the Poincare 
generators, 

P l =p\ M*j = x l pi - X ip i , ( . 

P°=H t , M M = x l H t . [W) 

This outcome is as expected: Compared to (54), p° has been replaced by H t , and 
a; has been set to zero. It should, however, be pointed out that for non-Cartesian 
coordinates the construction of the Poincare generators is less straightforward. 

Let us address the question of kinematical versus dynamical generators. In 
agreement with (58) and (59), one has 

N* = = x l N j - x j N i , i,j = 1,2,3, (61) 

so that S is both translationally and rotationally invariant confirming that the 
dimension of its stability group is six (cf. Table 1). On the other hand, 

iV° = 1 ^ , (62) 
x°N l - x l N° = -x l ^ , (63) 

from which we read off that, apart from the Hamiltonian, also the boosts are 
dynamical, i.e., S is not boost invariant. The latter fact is, of course, well known 
because the boosts mix space and time. Under a boost along the n-direction with 
velocity v, t transforms as 



t — > t' = t cosh w + (n • x) sinh u> , 



(64) 
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where n = v/v and w is the rapidity, denned through tanhw = v. From (64) it 
is evident that the hypersurface S : t = is not boost invariant. 

In obtaining the representation (60), wc make the Poincare algebra com- 
patible with the instant-form gauge-fixing constraint, x° = 0. An elementary 
calculation, using {x l , p J } = <5 y , indeed shows that the generators (60) really 
obey the bracket relations (52). We have already seen in (49) that the Hamilto- 
nian P° = H t generates the correct dynamics. 

At this point it is getting time to really consider an alternative to the instant 
form in some detail. 



2.3 The Front Form 

For an arbitrary four-vector a we perform the following transformation to light 
cone coordinates, 

(a , a 1 , a 2 , a 3 ) (a + , a 1 , a 2 , a~) , 



where we have defined 



a + a 3 



We also introduce the transverse vector part of a as 

a± = (a 1 , a 2 ) . 

The metric tensor (29) becomes 



h 



a0 



n 2 h 1 
h -H 



( 

-1 

Vl/2 



1/2\ 

-1 
0/ 



(65) 



(66) 



(67) 



(68) 



The entries 1/2 imply nonvanishing h and thus a slightly unusual scalar product, 
a- b = g^b" = \a + b- + \ a -b+ - a l b l , i = 1,2 . (69) 



According to Table 1, the front form is defined by choosing the hypersurface 
E : x + = 0, which is a plane tangent to the light-cone. It can equivalently be 
viewed as the wave front of a plane light wave traveling towards the positive 
^-direction. Therefore, S is also called a light-front. The normal vector is 



N= (1,0,0,-1) , N 2 = Q, 



(70) 



where TV has been written in ordinary coordinates. We see that N + = N° + 
N 3 = which implies that the normal N to the hypersurface lies within the 
hypersurface (Neville and Rohrlich 1971; Rohrlich 1971). As A is a light-like or 
null vector, S is often referred to as a null-plane (Neville and Rohrlich 1971; 
Coester 1992). We have depicted the front-form hypersurface £ together with 
the light-cone in Fig. 1. 
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Fig. 1. The hypersurface E : x + — defining the front form. It is a null-plane tangential 
to the light-cone, x 2 = 0. 




As is already obvious from (68), the unit vector in x + -direction is another 
null-vector, 

^=9^ = ^,0,1), (71) 

so that n ■ N = 1 as it should. Given the scalar product (69), we infer the 
invariant distance element 

ds 2 = g^dx^dx" = dx + dx~ - dx i dx i = ( ^ dxdx \ dx+dx+ 

\ dx + dx + dx + ) 

from which the einbein h can be read off as 

h(x+) = ±- - x L x l = v- - vi . (73) 

Note that velocities are dimensionless, so that despite appearance the result is 
consistent (if you do not like it as it stands, just insert the appropriate factors 
ofc). 
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The Hamiltonian is obtained by solving the constraint p + p — p\ — m 2 = 0, 
which is now linear in p~ . The result is 

H x+ =n.p = p-/2 = ^±^. (74) 

Let me reemphasize that this Hamiltonian does not contain a square root as 
already pointed out by Dirac. However, now it is crucial that the FP operator 
is nonvanishing, 

FP = —2N ■ p = -2p + 7^ . (75) 

While this is always true for massive particles, it is violated for massless 'left- 
movers', i.e. for particles travelling in the negative z-direction at the speed of 
light. In this case, we have a 'Gribov problem' (Gribov 1978), as the particles 
move within our gauge-fixing hyperplane, x + = 0. We will return to this issue 
later on. 

For massive particles, the dynamics is consistently generated by means of the 
Poisson brackets 

v- =x~ ={x~ ,H x+ } = ^ , v i =x i = {x i ,H x+ } = ^ . (76) 

Note finally, that the Hamiltonian (74) is not the normal projection N ■ p of the 
momentum, because N ■ p lies within S and thus corresponds to a kinematical 
direction. 

As for the instant form, the light cone representation of the Poincare gen- 
erators can be obtained by solving the constraint p 2 = m 2 for p~ , inserting 
the result into the elementary representation (54) of the generators and setting 
x + = 0. The kinematical generators are 

P*=p\ P + =p + , 
M +t = -x'p + , M 12 = x 1 p 2 - x 2 p 1 , M+- = -x-p + . (77) 

They correspond to transverse and longitudinal translations within S (P l and 
P + , respectively), transverse boosts and rotations (M +l ), rotations around the 
z-axis (M 12 ) and boosts (!) in the z-direction (M + ~). The latter will be further 
analysed in a moment. We thus have found seven kinematical generators, so 
that the front form leads to the largest stability group among Dirac's forms of 
dynamics (cf. Table 1). 

The dynamical generators are given by 

P~ = p2 t +m2 , M- 1 = x-p l - x l p- . (78) 
p+ 

As expected, the M~ l depend on the Hamiltonian, p~ . If we now consider rota- 
tions around the x- or y-axis, generated by 

L \ =M 23 = i(M 2 +-M 2 ") , (79) 

L 2 = M 31 = i(M +1 — M _1 ) , (80) 
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we note that they correspond to dynamical operations due to the appearance 
of M~ l . This leads to the notorious 'problem of angular momentum' within 
the front form, see e.g. (Fuda 1990). Except for the free theory, it is very hard 
to write down states with good angular momentum as diagonalizing L 2 is as 
difficult as solving the Schrodinger equation. A similar problem arises for parity. 
This exchanges light-cone space and time and thus also becomes dynamical 
(Burkardt 1996). For the kinematical component of the angular momentum, 
L z = M 12 , these difficulties do not arise. 

Consider now the following boost in z-direction with rapidity lo written in 
instant-form coordinates, 

t — > t cosh lo + z sinh lo , (81) 
z — > t sinh oj + z cosh lo . (82) 

As stated before, such a boost mixes space and time coordinates z and t. If we 
add and subtract these equations, we obtain the action of the boost for the front 
form, 



(83) 
(84) 



We thus find the important result that a boost in z-direction does not mix 
light-cone space and time but rather rescales the coordinates! Note that x + and 
x~ are rescaled inversely with respect to each other. The scaling factor can be 
written as 

- V^- (85) 

if the rapidity lo is defined in the usual manner in terms of the boost velocity v, 
tanhw = v. One should note in particular, that one has the fixed point hyper- 
surface S : x + — which is mapped onto itself, so that the relevant generator, 
M + ~ — 2M 30 = —2K 3 , is kinematical, in agreeement with (77). However, we 
see explicitly that this is no longer true for x + ^ 0, where we get a rescaling of 
x + . Stated differently, the transformation to light-cone coordinates diagonalizcs 
the boosts in ^-direction. Therefore, the behavior under such boosts becomes 
especially simple. A pedagogical discussion and some elementary applications 
can be found in (Parker and Schmieg 1970). 

We are actually more interested in the transformation properties of the mo- 
menta, as these, being Poincare generators, are more fundamental quantities 
than the coordinates, in particular in the quantum theory (Leutwyler and Stern 
1978). As transforms as a four-vector we just have to replace x^ by P M in 
the boost transformations (83, 84) and obtain, 



P+ -» e w P+ , 
P~ -» e~ u P- , 



(86) 
(87) 
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We remark that P + — is a fixed point under longitudinal boosts. In quantum 
held theory, it corresponds to the vacuum. For the transverse momentum, P l , 
one finds a transformation law reminiscent of a Galilei boost, 

P l -» P l + v l P + . (88) 

In this identity, describing the action of M +t , longitudinal and transverse mo- 
menta (which arc both kinematical) get mixed. 

We can now ask the question how to boost from (P + ,P 1 ) to momenta 
(Q + , Q l ). This can be done by hxing the boost parameters u) and v l as 

W = -Iogg, v* = 2^£-. (89) 

Obviously, this is only possible for P + ^ 0. We emphasize that in the construc- 
tion above there is no dynamics involved. For the quantum theory, this means 
that we can build states of arbitrary light-cone momenta with very little effort. 
All we have to do is applying some kinematical boost operators. The simple be- 
havior of light-cone momenta under boosts will be important for the discussion 
of bound states in Section 4. 

The similarity between (88) and Galilei boosts is not accidental. This is ex- 
hibited by the following subalgebra of the light-cone Poincare algebra (Susskind 
1968; Bardakci and Halpern 1968). Consider the Poisson bracket relations of the 
seven generators P^, M 12 , M + \ 

{M 12 , M +t } =e ij M +j , 
{M 12 , P 1 } = e ij P j , 
{M +l ,P~} = -2P i , 
{M +i , Pi) = -6*3 P+ 

All other brackets of these generators vanish. Compare now with the two- 
dimensional Galilei group. Its generators (for a free particle of mass (i) are: two 
momenta k % , one angular momentum L = e IJ f y , two Galilei boosts G l — (ix l , 
the Hamiltonian H = k l k l /2(i and the mass /i, which is the Casimir genera- 
tor. Upon using {x l ,k j } = 6 ij and identifying P l <-» k\ M 12 <-> L, M +i <-» 
— 2G\ P + <-» 2(i and P~ <-» H, one easily finds that (90) forms a subal- 
gebra of the Poincare algebra which is isomorphic to the Lie algebra of the 
two-dimensional Galilei group. (A second isomorphic subalgebra is obtained via 
identifying M~ % <-> 2G % and exchanging P + with P~ .) The first two identi- 
ties in (90), for instance, state that M +l and P l transform as ordinary two- 
dimensional vectors. P + can be interpreted as a variable Galilei mass which is 
also obvious from the nonrelativistic appearance of the light-cone Hamiltonian, 
P~ = (Pl + m 2 )/P+ and the Galilei boost (88). 

One thus expects that light-cone kinematics will partly show a nonrelativistic 
behavior which is associated with the transverse dimensions and governed by the 
two-dimensional Galilei group. This expectation is indeed realized and leads, 



(90) 
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for instance, to a separation of center-of-mass and relative dynamics in multi- 
particle systems. This will be discussed at length in the beginning of Section 4. 

So far, our discussion of the Poincare algebra was restricted to the free case. 
With the inclusion of interactions, one expects all dynamical Poincare generators 
to differ from their free counterpart by some 'potential' term W. This has already 
been pointed out by (Dirac 1949), who also stated that finding potentials which 
are consistent with the commutation relations of the Poincare algebra is the "real 
difficulty in the construction of a theory of a relativistic dynamical system" with 
a fixed number of particles. 

It has turned out, however, that Poincare invariance alone is not sufficient 
to guarantee a reasonable Hamiltonian formulation. There are no-go theorems 
both for the instant (Leutwyler 1965) and the front form (Jacn ct al. 1984), 
which state that the inclusion of any potential into the Poincare generators, 
even if consistent with the commutation relations, spoils relativistic covariance. 
The latter is a stronger requirement as it enforces particular transformation laws 
for the particle coordinates. Thus, covariance imposes rather severe restrictions 
on the dynamical system (Leutwyler and Stern 1978). 

The physical reason for these problems is that potentials imply an instan- 
taneous interaction-at-a-distance which is in conflict with the existence of a 
limiting velocity and retardation effects. Relativistic causality is thus violated. 
This is equivalently obvious from the fact that a fixed number of particles is in 
conflict with the necessity of particle creation and annihilation and the appear- 
ance of antiparticles. 

Nevertheless, with considerable effort, it is possible to construct dynamical 
quantum systems with a fixed number of constituents which are consistent with 
the requirements of Poincare invariance and relativistic covariance (Leutwyler 
and Stern 1978; Sokolov and Shatnii 1979; Coestcr and Polyzou 1982). 

At this point one might finally ask whether the different forms of relativistic 
dynamics are physically equivalent. From the point of view that different time 
choices correspond to different gauge fixings it is clear that equivalence must 
hold. After all, we are just dealing with different coordinate systems. People 
have tried to make this equivalence more explicit by working with coordinates 
which smoothly interpolate between the instant and the front form (Prokhvatilov 
and Franke 1988; Prokhvatilov and Franke 1989; Lenz et al. 1991; Hornbostcl 
1992). In the context of relativistic quantum mechanics, it has been shown that 
the Poincare generators for different forms are unitarily equivalent (Sokolov and 
Shatnii 1979). 

We are, however, more interested in what might be called a 'top-down ap- 
proach'. Our aim is to describe few-body systems not within quantum mechanics 
but quantum field theory to which we now turn. 
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3 Light-Cone Quantization of Fields 
3.1 Construction of the Poincare Generators 

We want to derive the representation of the Poincare generators within field 
theory and their dependence on the hypersurface S chosen to define the time 
evolution. To this end we follow (Fubini et al. 1973) and describe the hypersur- 
face mathematically through the equation 

S : F(x) = t . (91) 

The surface element on S is implicitly defined via 

J da tl u{x)= j d i xN^5{F{x)-T)u{x) , (92) 

where, as before, — d^F(x) is the normal on S and u some integrable 
function. We will write this expression symbolically as 

da^ = d A x N»5{F{x) - t) . (93) 

The central object of this subsection will be the energy-momentum tensor, 

dC 

9{9^4>) 

with C being the Lagrangian depending on fields that are collectively denoted by 
<f>. With the help of the energy-momentum tensor (94) we can define a generator 

A[f] = J^f v (x)T^(x) , (95) 

where A and / can be tensorial quantities carrying dummy indices a, (3,. . . which 
we have suppressed. A[f] generates the infinitesimal transformations 

5 f B(x)=f4x)d»B(x), (96) 
5 f x" - F(x) , (97) 

where / is now understood as being infinitesimal. In the same way as for a finite 
number of degrees of freedom, the generator A is called kinematical, if it leaves 
S invariant, that is, 

6 f F = f^F = f-N = 0. (98) 

Otherwise, A is dynamical. With the energy-momentum tensor at hand, 
we can easily show that kinematical generators are interaction independent. We 
decompose T Miy , 

= T» v - g^C int , (99) 
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into a free part Tq V = T^ v (g = 0), g denoting the coupling, and an interacting 
part (we exclude the case of derivative coupling). If A is kinematical, we have 
from (93, 95, 98), 

A int [f] = -J d A x5(F-T)C mt f»d^F = - J d A x5{F-T)C- mt 8 f F = 0, (100) 

which indeed shows that A does not depend on the interaction. Dynamical op- 
erators, on the other hand, will contain interaction dependent pieces. Of course, 
we are particularly interested in the Poincare generators, P a and M a/3 . They 
correspond to the choices f£ = g^ and f^ 13 — x a g^~x l3 g^, respectively, so that, 
from (95), they are given in terms of T^ v as 



da^ a , (101) 

M aP = J da^{x a T^ 3 - x' 3 T^ a ) . (102) 

From (98) it is easily seen that the Poincare generators defined in (101, 102) 
act on F(x) = r exactly as described in (57). The remarks of Section 2 on the 
kinematical or dynamical nature of the generators in the different forms are 
therefore equally valid in quantum field theory. 

Let us first discuss the instant form. We recall the hypersurface of equal time, 
S : F(x) = N ■ x = x° = t, which leads to a surface element 

do* =d A xN^8{x a -t) , JV" = (1,0). (103) 

Using (101, 102), the Poincare generators are obtained as 

P» = J d 3 xT ^ , (104) 

M"" = j d 3 x (x»T Qv - x u T a ») . (105) 

For the front form, quantization hypersurface and surface element are given by 

£ : F(x) = N ■ x = x + = r , da^ = d A x 5(x + - r) , (106) 

where N is the light-like four- vector of (70). In terms of T M1/ , the Poincare 
generators are 

P" = ± J dx-d 2 x ± T +f * , (107) 

= \\ dx-d 2 x ± (x»T + " - x u T + ») . (108) 

J E 

The somewhat peculiar factor 1/2 is the Jacobian which arises upon transforming 
to light-cone coordinates. 
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3.2 Schwinger's (Quantum) Action Principle 

Our next task is to actually quantize the fields on the hypersurfaces E : r = 
F(x) of equal time r. There is more than one possibility to do so, and we 
will explain a few of these. We begin with a method that is essentially due 
to Schwinger (Schwinger 1951; Schwinger 1953a; Schwinger 1953b). We define a 
four-momentum density 

dC 

so that the energy-momentum tensor can be written as 

T i*v = md v <t>-g» v C. (110) 

In some sense, this can be viewed as a covariant generalization of the usual Leg- 
cndre transformation between Hamiltonian and Lagrangian. Using the normal 
of the hypersurface E, we define the canonical momentum (density) as the 
projection of 27^, 

n = N-n. (Ill) 

Schwinger's action principle states that, upon variation, the action S = j d d xC 
changes at most by a surface term which (if E is not varied, i.e. 5x^ = 0) is 
given by 

SG{t) = j da„n»5cf) = J d d x5(T- F)nS<f> . (112) 

The quantity SG is interpreted as the generator of field transformations, so that 
we have 

6<p = {</>, SG} , (113) 

in case that E is entirely space-like (with time-like normal) (Schwinger 1951; 
Schwinger 1953a). We note in passing that the generator SG is a field theoretic 
generalization of the canonical one-form dG = Pidq 1 used in analytical mechan- 
ics. 

As in the preceding section we have to distinguish two cases depending on 
whether the normal vector N is time-like or space-like. For time-like N, the 
associated hypersurface is space-like. The basic example for this case is the 
instant form, to which we immediately specialize. The canonical momentum 
density is given by the velocity, tt — <f>, and the Lagrangian is quadratic in <fi. 
The canonical Poisson bracket is derived from Schwinger's action principle using 
(113), 

54>(x) = {<t>(x) , SG(t)} 

= jdy° j cPy8{x -y°){<f>(x),n{y)}8<l>(y) 

= J d 3 y {fa), 4>{y)}5<t>{y)\ x0=y0=T . (114) 
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The canonical Poisson bracket, therefore, must be 

{<t>(x),<f>(y)} x> yei; = 5 3 (x-y)> 

which, of course, is the standard result. The second case, N light-like, corre- 
sponds to the front form. With minor modifications, Schwingcr's approach can 
also be used here, resulting in what is interchangably called light-cone, light- 
front or null-plane quantization. The canonical light-cone momentum is 

7T = N ■ n = N ■ d<f> = d + (t> = 2-^—<j> , (116) 

ox~ 

which is peculiar to the extent that it does not involve a (light-cone) time deriva- 
tive. Therefore, 7r is a dependent quantity which does not provide additional 
information, being known on S when the field is known there. Thus, ir is merely 
an abbreviation for d + <f> which is a spatial derivative. Again, the reason is that 
the normal of the null-plane S lies within S. This leads to the important 
consequence that the light-cone Lagrangian is linear in the velocity <9~</>, or, put 
differently, that light-cone field theories are first-order systems. As a result, 4> 
and d + 4> have to be treated on the same footing within Schwingcr's approach 
which leads to an additional factor 1/2 compared to (113), 

±64> = {<(>, 6G} , (117) 

with a front-form generator 

5G{x+) = \J dx-d 2 x ± d + cf>5<f) . (118) 

The appearance of the peculiar factor 1/2 in (117) has been discussed at length 
by (Schwinger 1953b) — see also (Chang et al. 1973). Roughly speaking it stems 
from the fact that the independent field content within the front form is only one 
half of that in the instant form. The factor 1/2 cancels the light-cone Jacobian 
J = 1/2 in (118), so that we are left with the Poisson bracket, 

{4>(x) , n(y)} x yes = S(x~ - y-)S 2 ( Xl _ - yj . (119) 

As usual, commutators are inferred from Poisson brackets by invoking Dirac's 
correspondence principle, that is, by replacing the bracket by i times the com- 
mutator. For arbitrary classical observables, A, B, this means explicitly, 

[A,B]=i{aTb} . (120) 

We do not address the question of operator-ordering ambiguities at this point, 
as these will not be an issue in the applications to be discussed later on. One 
should, however, be aware of this problem, as it indeed can arise within the 
framework of light-cone quantization (Hcinzl et al. 1996). 

Using (120), the bracket (119) leads to the following commutator, 



[<j>{x), ir{y)] x+=y+=T = i8(x -y )5 2 (xj_ - y_J 



(121) 
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As the independent quantities are the fields themselves, we invert the derivative 
d + and obtain the more fundamental commutator 

[<Kx),<Kv)]x+=v+=t = -isgn(ar - 2T)<S 2 (x ± - yj . (122) 

In deriving (122) we have chosen the anti-symmetric Green function sgn(ar) 
satisfying 

^-sgn(x-) = 2d(x-), (123) 

so that (121) is reobtained upon differentiating (122) with respect to y~. We 
will see later that the field commutator (122) can be derived directly within 
Schwingcr's method. Before that, however, let us study the relation between the 
choice of initializing hypersurfaces, the problem of field quantization and the 
solutions of the dynamical equations. 



3.3 Quantization as an Initial and/or Boundary Value Problem 

As a prototype field theory we consider a massive scalar field (p in 1+1 dimen- 
sions. Its dynamics is encoded in the action 

S{(j)} = J d 2 xC = j d 2 x(^d l ,(f>d^<f)- \m 2 (t) 2 - V[<f>]) , (124) 

where V is some interaction term like e.g. A</> 4 and C = Co + V. By varying the 
free action in the standard way we obtain 



5S= f da t ,n"6<j)+ [ 

JdM JM 



(125) 



If we do not vary on the boundary of our integration region M, S(j)\gM = 0, the 
surface term in SS (which is closely related to 5G from (112)), vanishes and we 
end up with the (massive) Klein-Gordon equation in 1+1 dimensions, 

(□ + m 2 )(/> = 0. (126) 

In this subsection, we will solve this equation by specifying initial and/or bound- 
ary conditions for the scalar field <j) on different hypersurfaces S. In addition, 
we will clarify the relation between the associated initial value problems and the 
determination of 'equal-time' commutators. 

It may look rather trivial to consider just the free theory, but this is not 
entirely true. Let us analyze what quantization of a field theory means in the 
light of the different forms of relativistic dynamics. One specifies canonical com- 
mutators like [4>{x), <j>(y)]x,y£x:, where the hypersurface E: r = const defines the 
evolution parameter r. As both x and y lie in S, the commutator is evaluated 
at 'equal time', which implies that it is a kinematical quantity. Therefore, it is 
the same for the free and the interacting theory. 
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Now, if <j) is a free field, the commutator, 



[<f>(x),<f>(0)} = %A{x) , 



(127) 



is exactly known: it is the Pauli-Jordan or Schwinger function A (Jordan and 
Pauli 1928; Schwinger 1949) which is a special solution of the Klein-Gordon 
equation (126). It can be obtained directly from the action in a covariant manner 
as a Peierls bracket (Peierls 1952; DeWitt 1984). Alternatively, one can find it 
by evaluating the Fourier integral, 



where I have given both the instant and front form representation (Heinzl and 
Werner 1994). We note that A is antisymmetric, A(x) = —A(—x) and Lorentz 
invariant (under proper orthochronous transformations). Most important, it is 
causal, i.e. it vanishes outside the light-cone, x 2 < (see Fig. 2). 



Fig. 2. The Pauli-Jordan function as a function of T = mx /2 and X = mx 1 /2. It 
vanishes outside the light-cone and oscillates inside. 




±sgn(x°)6{x 2 ) J (mV^2) 



|[sgn(x + ) + sgn(x )] J (mVx+x-) , 



(128) 
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If (j) is an interacting field, causality, of course, must still hold. If x and y are 
space-like with respect to each other, the commutator thus still vanishes, 

[0(aO,0(i/)](*-„)'<o = O. (129) 

This expresses the fact that fields which are separated by a space-like distance 
cannot communicate with each other. For the front form, with the hypersurface 
S : x + = 0, (129) cannot be used to obtain the canonical commutators: In 1+1 
dimensions, S is part of the light-cone and therefore entirely light like. In higher 
dimensions, £ still contains light like directions namely where x~ = x± = 0. For 
this reason, the light front commutator (122) of two free fields does not vanish 
identically. 

Let me now discuss the explicit relation between the choice of equal-time 
commutators and the classical initial/boundary for the Klein-Gordon equation. 
Three examples are of interest. 



Cauchy Data: Instant Form. Conventional quantization on a space-like sur- 
face (based on the instant form) corresponds to a Cauchy problem: if one specifies 
the field </> and its time derivative <j) on £ : x° = t = 0, 

0(t = O, *) = /(*), (130) 
j>(t = 0,x)=g(x) , (131) 

where the functions / and g denote the initial data (depending on x = x 1 ), 
the solution of the Klein-Gordon equation is uniquely determined. This can be 
checked by considering the Taylor expansion around (i, x) € £, i.e. t — 0, 

cp(t, x) = 0(0, x) + i0(O, x) + ±t 2 0(O, x) + . . . , (132) 

with the overdot denoting the time derivative. From this we see that one has to 
know all time derivatives of <f> on E once the data /, g are given. If we calculate 
these, 

<£ = /, 

ij> = 4>"- m 2 ^ = f" - m 2 f , 

^ 3< ^ ill 2 2 II 2 

— = -m 4> = g -m g , 

i . (133) 

we find that indeed all time derivatives are given in terms of / and g and their 
known spatial derivatives, denoted by the prime. In the last two identities, we 
have made use of the equation of motion. As a result, we see that the Cauchy 
problem is well posed: the solution of the Klein-Gordon equation is uniquely 
determined by the data on S. 
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Upon quantization, this translates into the fact that the Fock operators can 
be expressed in terms of the data, 

aip 1 ) = J dx 1 e- iplxl u p (j)(x = 0,x 1 )+i<p{x° = 0,x r ) , (134) 

with uj p — (jp\ + m 2 ) 1 / 2 . In addition, the canonical commutators can be viewed 
as the Cauchy data for the Pauli- Jordan function A, 



4>{x),<f>(0) 

<P(x),<j>(0) 



x°=0 



iA(x)\ x0=0 = , 
= iM x )\x°=o = -iSix 1 ) . 



(135) 
(136) 



As stated above, the vanishing of the commutator (135) is due to causality. 



The Characteristic Initial- Value Problem. In the following I will perform 
an analogous discussion for the hypersurfaces £ : x^ — const, which, in d = 1+1, 
constitute the entire light-cone, x 2 = 0. In Dirac's classification, the light-cone 
corresponds to a degenerate point form with parameter a = (see Table 1). 
One thus does not have transitivity as points on different 'legs' of the cone are 
not related by a kinematical operation. Still, it turns out that the associated 
initial value problem is well posed (Neville and Rohrlich 1971; Rohrlich 1971) 

The light-fronts x^ = are characteristics of the Klein-Gordon equation 
(Domokos 1972). Therefore, one is dealing with a characteristic initial- value 
problem (Courant and Hilbert 1962), for which one has to provide the data 



<j>(x + =0,x~) = f(x-) , 
4>{x + ,x- =0) =g(x + ) , 
f(x~ = 0) = g(x + = 0) 



(137) 
(138) 
(139) 



where the last identity is a continuity condition. Consistency is again checked 
by Taylor expanding, this time around (x + , x~) = 0, 



d+4 


> = d + .f = 


/', 




> = d~g = 


9 , 




> = m 2 <p — 


m 2 f 




> = /", 






' = 9 , 






> = m 2 f, 




d'd-d+4 


) = m 2 g , 





m 2 g 



; . (wo) 

Whereever a factor of m 2 appears we have made use of the Klein-Gordon equa- 
tion. We thus note that the data (together with their known derivatives) de- 
termine all partial derivatives of <f> at the vertex of the cone, (x + ,x~) = 0. 
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Intuitively, this corresponds to the fact that the information spreads from a 
source located at origin. 

The characteristic initial- value problem amounts to quantization on two char- 
acteristics, x ± = 0, i.e., in d = 1 + 1, really on the light cone, x 2 = 0. The 
following two independent commutators, 

[4>{x), 0(O)] I±=O = t4(aOU ±=0 = - jsgn(arF) , (141) 

are then characteristic data for the Pauli- Jordan function. It turns out that, in 
case the field <p is massless, the above quantization procedure is the only con- 
sistent one (in d=l+l), if one wants to use light-like hypersurfaces (Bogolubov 
et al. 1990). 

However, it is important to note that the characteristic initial- value problem 
does not correspond to light-cone quantization. One would need two Hamilto- 
nians P~ and P + , and, accordingly, two time parameters. This seems somewhat 
weird, to say the least, and will not be pursued any further. 



Initial-Boundary-Data: Front-Form. In order to find the initial-value prob- 
lem of the front form with a single time parameter x + , let us naively try a 
straightforward analog of the Cauchy data and prescribe field and velocity on 
E : x+ = 0, 

4>{x + =0,x-) = f(x-), (142) 
d-(j>{x + = 0,aT) =g(x~) . (143) 

It turns out, however, that this overdetermines the system. Namely, from the 
equation of motion 

d + d-(f> = m 2 (t> , (144) 

it is actually possible to obtain the velocity d~<p by inversion of the spatial 
derivative <9 + , 

d-(j> = m 2 {d + )- 1 4> = m 2 {d + y 1 f . (145) 

The last identity holds on U and implies that the data g are unnecessary (and 
will even lead to an inconsistency) as the velocity is already determined by /. 
This is confirmed by the Taylor expansion on E, 

<i> = f, 
= d+f , 

d-4 - m 2 (d + )-\f , 
d~d-(f> = m 2 (d+)- 1 d-(f> = m 4 (<9+)- 2 / , 

; . (U6) 

It thus seems that the front form requires only half of the data as compared to 
the instant form. This appearance, however, is deceptive. Note that we have to 
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invert the differential operator d + . The inverse is nothing but the Green function 
G defined via 

d + G(x~) = 6{x~) . (147) 

Clearly, this Green function is determined only up to a homogeneous solution h 
satisfying 

d+h = , (148) 

i.e. up to a zero mode h = h{x + ) of the operator <9 + . Thus, in order to uniquely 
specify the Green function (147), we have to provide additional information in 
terms of boundary conditions. The standard choice is to demand antisymmetry 
in i", whence 

G(x-) = hgn(x-) , (149) 

which we have already used in (122) and (123). Before we discuss the physi- 
cal reason for demanding antisymmetry, let us briefly go to momentum space 
where we replace d + by ip + . The equation (147) for the Green function becomes 
ip + G(p + ) — 1, which has the general solution 

G(p + ) = -i/p+ + h{p-)6(p+) . (150) 

In this identity, l/p + has to be viewed as a distribution corresponding to an 
arbitrary regularization of the singular function l/p + (Gelfand and Shilov 1964). 
Any two regularizations differ by terms proportional to S(p + ), i.e. a zero mode of 
p + . Choosing an antisymmetric Green function uniquely yields a principal value 
prescription, 

iG{p + ) = V— = \( — - — + — - — , (151) 
KL ' p+ 2\p++iep + -ie)' K ' 

which is the canonical regularization of l/p + . 

Altogether we have seen that the front form corresponds to prescribing both 
initial and boundary conditions, so that one has a 'mixed' or initial-boundary 
value problem. What are the implications for quantization? We address this ques- 
tion by determining the Poisson brackets through the requirement that Euler- 
Lagrange and canonical equations should be equivalent. To this end we solve 
the Klein-Gordon equation (126) for the velocity 4> = d(j)/dx + as in (145). This 
gives 

777 / 

ct>(x+,x-) = -— J dy-G{x-,y-)<j>{x+,y-) , (152) 

using the Green function G from (149). The Hamiltonian equation of motion is 
given by the Poisson bracket with H = \m 2 J dx~(f> 2 , 

m 2 f 

<f>(x+,x-) = — I dy~ {4>(x + ,x-),<f>(x+,y-)}<f>(x+,y-) (153) 

with the bracket of the fields 4> to be determined. Clearly, Euler-Lagrange and 
Hamiltonian equation of motion, (152) and (153) become equivalent if one iden- 
tifies 

{<j>{x + ,x~), (j>{x + ,y~)} = -\G{x~ ,y~) . (154) 
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We thus see that the fundamental Poisson bracket coincides with the Green 
function which, accordingly, justifies the requirement of antisymmetry. After 
quantization, (154) of course coincides with (122), the result from Schwinger's 
action principle, specialized to d = 1 + 1, 

[<t>{x), <f>(0)] x += = iA{x)\x+=o = -^ s g n 0O • (155) 
From the momentum space perspective, 

W + ),m\ = {G{p+)= l -V^, (156) 

we conclude that, technically, light-cone quantization is the inversion of the lon- 
gitudinal momentum p + and as such requires the specification of initial-boundary 
data. In some sense, this can also be viewed as an infrared regularization because 
one provides a prescription of dealing with a pole at vanishing longitudinal mo- 
mentum, p + — 0. As is well known, a particularly nice way of regularizing in the 
infrared is to enclose the system under consideration in a finite spatial volume. 
This is the topic of the next subsection. 



3.4 DLCQ — Basics 

DLCQ is the acronym for 'discrctizcd light-cone quantization', originally devel- 
oped by (Maskawa and Yamawaki 1976; Pauli and Brodsky 1985a; Pauli and 
Brodsky 1985b) 4 . The physical system under consideration is enclosed in a fi- 
nite volume with discrete momenta and prescribed boundary conditions in x~ . 
Recently, there has been renewed interest in this method in the context of string 
theory (Banks et al. 1997; Susskind 1997; Lunin and Pinsky 1999). 

Our starting point is the Fourier representation for a solution of the Klein 
Gordon equation (still in infinite volume), 

4>{x)= J ^-X(p)6(p + p--m 2 )e-^ 

= j d ^ L x{p\p~)^{p- -m 2 /p + )e-^ 

f°° dn+ 

= / t^- [x(p + ,p-)e-^* + X (-P + , -p-y p - x ] 
Jo 47I 7? 

= / i^*(p + ) [ a (p + ) e ~** + a *(p + y p ' x ] ■ ( i57 ) 

The following remarks are in order: we have defined the on-shell energy, p~ = 
m 2 /p + ; contrary to the instant form, the integration over the positive and neg- 
ative mass hyperboloid is achieved by a single delta function. Again, this is a 

4 For a nice overview of recent developments see (Hiller 2000) . 
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consequence of the linearity of the mass-shell constraint in p~ . The two branches 
of the mass-shell correspond to positive and negative values of p + (and also p~ ) , 
respectively. Associated with the two signs of the kinematical momentum p + are 
the positive and negative frequency modes a, a* , defined in such a way that their 
argument p + is always positive (cf. the step function 9 in the last line) . This can 
cquivalcntly be viewed as the reality condition 

a*(p+) = a(-p+) , (158) 

as is obvious from the last step in the derivation (157). Upon quantization this 
implies that annihilation operators with negative longitudinal momentum p + are 
actually creation operators for particles with positive p + . The field commutator 
(155) is reproduced by demanding 

[a(k+), a\p + )} = 47rp+(5(fc + - p+) . (159) 

As already indicated, DLCQ amounts to compactifying the spatial light-cone 
coordinate, — L < x~ < i, and imposing periodic boundary conditions for the 
fields, 

(f)(x + ,x~ = —L) = </>(x + ,ar = L) , (160) 

which arc to hold for all light cone times x + . Space time is thus endowed with 
the topology of a cylinder. This implies discrete longitudinal momenta, fc+ = 
2wn/L, so that the Fock expansion (157) becomes 

<f>(x+ =0,x-)=a + y -L= (a n e- m ™~' L + a* n e in ™~/ L ) , (161) 

Note that we have allowed for a zero momentum mode ao- We will see in a 
moment that it actually vanishes in the free theory. Plugging (161) into the free 
Lagrangian 

L [<f>] = ^Jdx- (Ia+09-0- I m 2 2 ) , (162) 
we obtain (discarding a total time derivative) 

L [a n , a ] =-iy^ a n a* n - m 2 La\ - ~^~ a n a " = ~* ~ H > ( 163 ) 

n>0 n>0 n>0 

with H denoting the Hamiltonian and a* = da*Jdx + . From both representa- 
tions (162) and (163) it is obvious that the light-cone Lagrangian is linear in 
the velocity (d~<p and d* , respectively). A particularly suited method for quan- 
tization in this case is the one of Faddeev and Jackiw for first order systems 
(Faddeev and Jackiw 1988; Jackiw 1993). It avoids many of the technicalities 
of the Dirac-Bcrgmann formalism and is in general more economic. It reduces 
phase-space right from the beginning as there are no 'primary constraints' in- 
troduced. The method is essentially equivalent to Schwinger's action principle, 
especially in the form presented in (Schwinger 1953b). For the case at hand, the 
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method basically boils down to demanding equivalence of the Euler-Lagrange 
and Hamiltonian equations of motion (cf. last subsection) . 
The former are given by 

171?' L 

-ia n + a n = , (164) 

47m 

2m 2 La Q = . (165) 

The first equation, (164), is just the free Klein-Gordon equation which can be 
easily seen upon multiplying by fc+. The second identity, (165), is nondynamical 
and thus a constraint which states the absence of a zero mode for free fields, 
a = 0. 

The canonical equations are 

On = {a n , H} = ^2 ^-^ {a n , a* k } a k , (166) 

fc>0 

which obviously coincides with (164) if the canonical bracket is 

{a k , a* } = -iS kn ■ (167) 
The constraint (165) is obtained by differentiating the Hamiltonian, 

— = 2m 2 La = . (168) 
oa 

Let us briefly show that the approach presented above is equivalent to Schwinger's 
(Schwinger 1953b). From (163) we read off a generator 

5G = -iJ2a n Sa* n (169) 

n>0 

effecting the transformation 

da* n = {< , SG} = ~iJ2 K . *k} Sa* k , (170) 

fe>0 

which in turn implies the canonical bracket (167). 

Quantization is performed as usual by employing the correspondence princi- 
ple (120), so that, from (167), the elementary commutator is given by 

[a m ,al] = S mn . (171) 

The Fock space expansion for the (free) scalar field <f> thus becomes 

(f>(x + = 0, x~) = V —L= (a n e- m ™~' L + a n e in ™~/ L ) , (172) 
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Like in the infinite volume expression (157), the Fock 'measure' 1/V^Lm does 
not involve any scale like the mass m or the volume L. This is at variance with 
the analogous expansion in the instant form which reads 

*(m = o) = ^£-^^ . <™> 

where — L < x < L, k Obviously, the 'measure' 

(A; 2 + to 2 ) 1 / 2 does depend on m and L. We will discuss some consequences of 
this difference in Subsection 3.6. 

We can use the results (171) and (172) to calculate the free field commutator 
at equal light-cone time x + , 



(174) 



This coincides with (155) up to a finite size correction given by the additional 
term x~ /2L. The effect of this term is two-fold. First, it makes the sign function 
periodic (in the interval — L < x~ < L), and second, it guarantees the absence 
of a zero mode which must hold according to (165), (168), and 



I 



dx~[<l>(x),<l>(0)] x+=o = 0. (175) 



One may equally think of this as the finite-volume analog of the principal value 
prescription. 

The commutator (174) has originally been obtained in (Maskawa and Ya- 
mawaki 1976) using the Dirac-Bcrgmann algorithm for constrained systems. The 
Faddeev-Jackiw method, however, is much more economic and transparent. In 
particular, it makes clear that the basic canonical variables of a light-cone field 
theory are the Fock operators or their classical counterparts. The a n with, say, 
—N < n < N in (161) can be viewed as defining a (2N + l)-dimcnsional phase 
space. A phase space, however, should have even dimension. This is accomplished 
by choosing a polarization in terms of positions and momenta, here a n and a\, 
with n > 0, and by the vanishing of the zero mode, ao = 0. It turns out that 
this vanishing is a peculiarity of the free theory as is discussed in (Maskawa and 
Yamawaki 1976; Wittman 1989; Heinzl et al. 1992). 

At this point one should honestly state that the issue of zero modes is one 
of the unsolved problems of light-cone quantization 5 . The constraint equations 
for the zero modes arc in general very hard to solve unless one has some small 
parameter like in perturbation theory (McCartor and Robertson 1992; Heinzl 
et al. 1996) or within a large- expansion (Borderies et al. 1993; Borderies et al. 
1995). Using a path integral approach, it has recently been shown (Hellerman 

5 For an overview see e.g. Ch. 7 of (Brodsky, Pauli, and Pinsky 1998) and references 
therein. 
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and Polchinski 1999) that integrating out the zero modes constitutes a strong 
coupling problem. There are speculations that this problem might be less severe 
if one goes beyond quantum field theory, i.e. in string or M-theory (Banks 1999). 

In the last reference, the author also states that compactification in a light- 
like direction "is close to a space with periodic time" and thus "weird" , in view of 
possible 'grandfather paradoxes'. Therefore, the natural question arises whether 
DLCQ is actually consistent with causality. 

3.5 DLCQ Causality 

In this subsection I will address the question under which circumstances com- 
pactification of 'space' is compatible with the requirements of causality. The 
presented results are based on recent work with N. Scheu and H. Kroger (Heinzl 
et al. 1999). 

In (128) and (129) we have seen that the (infinite volume) commutator of 
two scalar fields vanishes whenever their space-time arguments are separated by 
a space-like distance (cf. Fig. 2). As already mentioned, this is a manifestation 
of the principle of microcausality, which is the general statement that the com- 
mutator of any two observables 0\(x) and 02 (j/) must vanish whenever their 
separation x — y is space-like. Physically, this implies that measurements of the 
observables 0\ and O2 performed at x and y, do not interfere. Some conse- 
quences of this principle are the spin-statistics theorem, analyticity properties 
of Green functions leading to dispersion relations etc. (Streater and Wightman 



Our starting point are the Fourier representations of the Pauli-Jordan func- 
tion, both for the instant and front form (denoted IF and FF, respectively), 



Both integrals yield the same result (128) for the Pauli-Jordan function. Note, 
however, that the integrand I(k + ) is exploding and rapidly oscillating for fc+ — > 
(see Fig. 3) so that the finite result for the integral is due to sizable cancellations 
that occur upon integration. 

To obtain the finite-volume representations for the commutator, one proceeds 
as in Subsection 3.4 by restricting the spatial coordinates, — L < x x ,x~ < L, 
and imposing periodic boundary conditions for the field (f>. Momenta become 
discrete, k\ = nn/L, and fc+ = 2ttti/L. The finite-volume representations are 
defined by replacing the integrals (176) and (177) by the discrete sums, 



1963). 




(176) 



(177) 




(178) 
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40 - 

20- 








0.3 


0.4 


k*/m 


-20 










-40 












N 1 

ra=l 


sin(fc„ • x) . 







The on-shell energies for discrete momenta are denned as 

uj n = (nV/L 2 + m 2 ) 1/2 and k~ = m 2 L/27m . (180) 

For both functions, Ajp and App, the periodicity in x 1 and x~ , respectively, 
with periodicity length 2L, is obvious. The limit N — > oo is understood unless 
we perform numerical calculations where N is kept finite. 

The evaluation of the sums (178) and (179) is not straightforward. To gain 
some intuition, we evaluate them numerically beginning with the IF expression 
(178). The resulting A IF is plotted in Fig. 4. 

Upon inspection, one notes the following: Up to small oscillations stemming 
from the (unavoidable) Gibbs phenomenon, Ajp vanishes outside the light- 
cone (|a;°| < < L), and thus is causal even in finite volume. If we let the 
summation cutoff N go to infinity, Ajp approaches the continuum Pauli- Jordan 
function A (for — L < a; , a; 1 < L). There is a clear physical picture behind 
these observations. One can imagine a periodic array of sources located at the 
quantization hypersurface x° = at points x 1 = 2Ln. These sources 'emit' 
spherical 'waves' into their own future LCs which start to overlap after time 
a; > L. At this point the 'waves' emanating from the sources begin to interfere. 
Thus, the influence of the BC is felt only after a long time (as large as the 
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Fig. 4. A IF (X,T) as a function of X = x 1 /2L. T = X° /2L = 0.2, mL = 1, N = 50. 
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spatial extension L of the system). This picture can be confirmed analytically. 
An application of the Poisson resummation formula yields 

A IF (x) =^2A(x°,x 1 +2Ln) , (181) 

n 

i.e. a periodic array of A's which are nonoverlapping as long as x° < L. 

For the front form, the situation turns out to be more complicated. Using 
Poisson resummation one can derive the finite-volume version of the canonical 
light cone commutator, at x + — 0, which is a periodic sign function, 

A FF (x + =0,x-) = sgn(x- + 2Ln) + x~/4L . (182) 

n 

This coincides with (174) if x~ is restricted to lie between — L and L. 

For x + 7^ 0, I have evaluated Aff numerically. The result is shown in Fig. 5 
as a function of the dimcnsionless variables v = x~ /2L, w = m 2 Lx + /2. For 
large values of w, Aff attains a very irregular shape, though numerically the 
representation (179) converges to a periodic function. The most important ob- 
servation, however, is that Aff does not vanish outside the light-cone, i.e. for 
x~ < 0, if x + > as in Fig. 5. This a clear violation of microcausality as has 
first been observed in (Scheu 1997). 

As already stated, it is not straightforward to confirm these findings analyti- 
cally. Poisson resummation does not work; first, because of the weak localization 
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Fig. 5. Aff(v,w) as a function of v = x /2L. w — 10000, N — 70. It does not vanish 
outside the LC, — 1 < v < 0. 
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properties of A in x~ (asymptotically, is goes like (x - )" 1 / 4 ); second, and even 
worse, because the zero mode I(k + — 0) does not exist for x + ^ 0. Nevertheless, 
an independent confirmation of causality violation can be obtained from resum- 
ming (179) in terms of Bernoulli polynomials, thereby replacing the Fourier series 
by a (rapidly converging) power series in w. The result is shown in Fig. 6 for 
w = 5. There is nice agreement with the Fourier representation (179) (and no 
Gibbs phenomenon, as expected). Again, causality violation is obvious. 

From a technical point of view, the violation of causality is not really aston- 
ishing. We are replacing the integral over the severely oscillating function I(k + ) 
by a Riemann sum with equidistant grid points. In this way, we are sampling the 
integrand in such a way that the huge cancellations present in the integral do 
not take place. Instead, for small k + , we replace I(k + ) by a random 'staircase' 
function which in the end produces the 'noise' seen in Fig. 5. 

At this point a natural questions arises: is there a remedy for the causality 
violation? The answer is positive. We have found two ways around the problem 6 , 
both, however, with shortcomings of their own. The first way is to regularize the 
integral (177), replacing A by A e in such a way that the associated integrand 
satisfies I e (k + — 0) = 0. One can chose e.g. a principal value regularization 
(Hornbostel et al. 1988) or a more sophisticated prescription (Nakanishi and Ya- 
mawaki 1977). In this way, one suppresses the oscillations and the divergence at 

6 For a third method, see (Chakrabarti et al. 1999). 
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Fig. 6. Comparison of the Fourier representation (179) with the result of Bernoulli 
resummation (smooth, heavy line). 




fc + = at the price of introducing a small causality violation of order e. But now 
A f can be approximated by a discrete sum if the momentum grid is sufficiently 
fine, Afc + -C e. The order of limits, however, becomes important. First, one has 
to perform the infinite volume limit, Afc+ — ► 0, L — > oo, and only then the 
limit e — > 0. For this method, Poisson resummation should work (Salmons et al. 
1999). However, it seems somewhat awkward and not very economic to perform 
two regularizations (finite L and e). 

An alternative way of resolving the problem is the following: instead of an 
equally spaced grid a la DLCQ (i.e. Afc+ = const) one can chose an adapted 
momentum grid with spacing A/c+ <~ 1/n for small n. In this way, one is sampling 
the small-fc + -region of I(k + ) in a more reasonable way. Practically, the method 
amounts to viewing Ajp as the correct finite-volume expression and replacing 
x° and x 1 by (x + ± x + )/2, respectively. This is equivalent to introducing new 
discrete momenta, = w n ± k\ . 

As a result, the point k + — becomes an accumulation point of the momen- 
tum grid which leads to a causal finite-volume representation A c (x + 7 x~) of A 
(see Fig. 7). This function, however, is no longer periodic in x~ . We thus find 
that, with a light-like direction being compactified, one cannot have both, peri- 
odicity and causality. On the other hand, the regularization method above seems 
to suggest that the causality violation is in some sense 'small' and thus may have 
a minor effect on the calculation of observables. Whether this statement is true 
has still to be worked out in detail. 
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Fig. 7. The causal commutator A c as a function of v. x + /2L = 0.2, mL = 50, N = 50. 




3.6 The Functional Schrodinger Picture 

Having discussed some difficulties of light-cone quantization in a finite 'box' we 
go back to the infinite volume and consider yet another method of quantization, 
the functional Schrodinger picture. The idea of this method is to mimic quantum 
mechanics within quantum field theory. States are described as functional &[<f>] 
depending on the field(s) <j>, while operators are combinations of multiplication 
by functional of <j> and functional differentiation with respect to <f). 



Bosons (Instant Form). For a free massive scalar field in two dimensions 
one finds the following results using instant form dynamics (Jackiw 1987). The 
canonical momentum (operator) acts via differentiation, 



8<j>{x) 



m ■ 



(183) 



The ground state $b of the system, the Fock vacuum, is obtained by direct anal- 
ogy with the harmonic oscillator. One rewrites the vacuum annihilation condition 
a(k)\0) — as a functional differential equation, 



/ 



dyf2(x - y)<j>{y) 



6(f>{x) 



M4>\ = o 



(184) 
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which is solved by a Gaussian, 

*bM~exp[-i(^/fy)] . (185) 
Here we have introduced the 'quadratic form', 

= J dxdycj>(x)n(x-y)<j>(y) , (186) 

using the kernel (or covariance) Q{x — y). The latter is defined by its Fourier 
transform, 

Q{k) = \/k 2 + m 2 = uo kl (187) 

which is nothing but the on-shell energy of a free massive scalar. For later use, 
it is important to note that the instant form covariance Q is explicitly mass 
dependent. As a consequence, if we have two free scalars of masses mi and mi, 
respectively, their Fock vacua are related by a Bogolubov transformation, 

|«2> = M^i) , (188) 
with the unitary operator U21 explicitly given by 

U21 = exp y" £ e k [ ai (k) ai {-k) - a\(k)a\(-k)} . (189) 

9k is the Bogolubov angle. We mention in passing that in order to properly define 
U21 as an operator one should work in finite volume to avoid infrared singularities 
(Heinzl 1998). Otherwise the two vacuum states have vanishing overlap. Within 
the functional Schrodinger picture this has been analysed by (Jackiw 1987). 



Fermions (Instant Form). For fermions, the situation is slightly more in- 
volved. Note that even within the instant form, the Dirac Lagrangian represents 
a first-order system, so that one expects some similarities with light -cone quan- 
tization in this case. This expectation will indeed turn out to be true. The instant 
form fcrmionic field operators are given by (Kiefer and Wipf 1994) 

** = 75(*" + 4)' ^ = 7i( < + £)' (190) 

and thus are linear combinations of multiplication by and differentiation with 
respect to the complex-valued Grassmann functions u a (x) and u* a {x). These 
functions characterize the states, for example the ground state (Fock vacuum) 
which is again Gaussian, 

&o[u,u*} ~ exp(w*,J?w) . (191) 
For 2d massive fermions, the covariance is found to be 

Q{k) = - 1 (fcg 1 - ma 3 ) , 
Vfc 2 + m 2 



(192) 
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with a 1 and cr 3 the standard Pauli matrices. Again, fi is explicitly mass depen- 
dent. In the massless case, m = 0, it becomes particularly simple, 

n(k) = sgn(fc)cr 1 , (193) 

or, after Fourier transformation, 

!2{x - y) = -V—^—a l . (194) 
7r x — y 

Here we have once more made use of the fact that the principle value P(l/x) is 
the Fourier transform of the sign function. 



Bosons (Front Form). As stated above, the latter case is somewhat similar to 
the generic situation in light-cone quantization. Let us again consider a massive 
free scalar field <j) m 2d with Fock expansion (157). We decompose it into positive 
and negative frequency part, 

(/) = (t) + [a] + (j)-[a*] = u + u* , (195) 

where <j)~~ = (4> + )* as <f> is real. Quantization is performed by defining field 
operators such that the canonical light-cone commutator (122) is reproduced. 
The solution turns out to be somewhat more complicated than for instant form 
fields, namely 

1 f S 

$+(x-) = u(x-) + -J dy-iA + {x- - V-)g^—j . (196) 

1 f s 

4>-(x-) - u*(x-) + -J dy-iA_(x- - y-)-^J—j ■ (197) 

A + and Z\_ are distributions that sum up to A, A + +A- = A, and are explicitly 
given by 

iA±(x~) = Tln(±a;- - ie) = tM x ~\ ± ^°(Tx~) ■ (198) 

The ground state is annihilated by <f) + , 4> +x l/ [u, u*] = 0, which yields a functional 
differential equation, again with Gaussian solution, 

i'oK u*] = cx P [-(m*, Qu)] . (199) 

The covariance is given by 

Q(x~) = 2id+6(x-) , i7(fc+) = 2fc+, (200) 

and thus is a local expression which is a very peculiar finding. In momentum 
space, the ubiquitous longitudinal momentum k + appears. One thing that is 
particularly obvious from (200) is the fact that the light-cone Fock vacuum 
is mass independent, &o(mi) — >^o( m 2) which means that the analog of the 
Bogolubov transformation (189) is trivial, i.e. U21 = 1. This has been checked 
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explicitly for several examples, including the Nambu-Jona-Lasinio model (Diet- 
maier et al. 1989) and bosons and fermions coupled to external sources (Heinzl 
1998). 

I believe that the locality of the covariance has far reaching consequences 
which, however, are still to be worked out in the present framework. Within 
ordinary Fock space language, some properties of the light-cone vacuum will be 
discussed in what follows. 

3.7 The Light Cone Vacuum 

One of the basic axioms of quantum field theory states that the spectrum of the 
four momentum operator is contained within the closure of the forward light- 
cone (Streater and Wightman 1963; Bogolubov ct al. 1975). The four-momentum 
of any physical, that is, observable particle thus obeys 

P 2 > , P° > , (201) 

which is, of course, consistent with the mass-shell constraint, p 2 = m 2 . The 
tip of the cone, the point P 2 = P° = 0, corresponds to the vacuum. From the 
spectrum condition (201) we infer that 

Po - Pi > Pi > or P° > \P 3 \ . (202) 

This implies for the longitudinal light-cone momentum, 

P+ = p° + p 3 > |p 3 | + p 3 > . (203) 

We thus have the important kinematical constraint that physical states must 
have nonnegative longitudinal momentum, 

(phys|P+|phys) > . (204) 

The spectrum of P + is thus bounded from below. Due to Lorentz invariance, the 
vacuum |0) must have vanishing four- momentum, and in particular 

P+|0) = . (205) 

Therefore, the vacuum is an eigenstate of P + with the lowest possible eigenvalue, 
namely zero. We will be interested in the phenomenon of spontaneous symmetry 
breaking, i.e. in the question whether — roughly speaking — the vacuum is degen- 
erate. Let us thus analyse whether there is another state, \p + = 0), having the 
same eigenvalue, p + — 0, as the vacuum. If so, it must be possible to create this 
state from the vacuum with some operator U, 

\p + - 0) = U\0) , (206) 

where U must not produce any longitudinal momentum. Note that within or- 
dinary quantization such a construction is straightforward and quite common, 
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for example in BCS theory. A state with vanishing three-momentum can be 
obtained via 

|p = 0) = J d 3 fc/(k) fl t( k ) a t(-k)|0) , (207) 

where / is an arbitrary wave function. Evidently, the contributions from modes 
with positive and negative momenta cancel each other. It is obvious as well, 
that within light-cone quantization things must be different as there cannot 
be an analogous cancellation for the longitudinal momenta which are always 
nonnegative. Instead, one could imagine something like 

|p+= Pi = 0) = J dk+J d 2 fc ± /(k ± )<5(fc + )a t (fc+,k ± )a t (fc + ,-k ± )|0) . (208) 

The problem thus boils down to the question whether there are Fock operators 
carrying light-cone momentum k + = 0. As we have seen in Subsection 3.4, there 
are no such operators, and a construction like (208) is impossible. 

The only remaining possibility is that, if U contains a creation operator 
a)(k + > 0) carrying longitudinal momentum k + ^ 0, there must be annihilators 
that annihilate exactly the same amount k + of momentum. Thus, after Wick 
ordering, U must have the general form 

U=(0\U\0)+ J dfc+/ 2 (fc + )a t (fc + )a(fc+) 

fc+>0 

+ J dp + J dk + f 3 (k + ,p + )a1( P + + k+)a(p + )a(k+) 

p+>0 fe+>0 

+ J dp + j dk + / 3 (fc + ,p + )a t (p + )a t (fc + )a(fc+ + p + ) 

p+>0 fc+>0 

+ . . . . (209) 

It follows that the light-cone vacuum |0) is an eigenstate of U, 

U\0) = (0|f/|0)|0) . (210) 

As we only deal with rays in Hilbert space, the action of U on the vacuum 
does not create a state distinct from the vacuum. One says that the vacuum 
is trivial. Put differently, 'there is no vacuum but the Fock vacuum'. Note that 
this is actually consistent with our findings in the last subsection: the light-cone 
vacuum is the same irrespective of the masses of the particles; the Bogolubov 
transformation present in the instant form becomes trivial. 

Let us analyse the dynamical implications of the general result (210). Any 
quantity that is obtained by integrating some functional of the fields over x~ , 
i.e., 

F[4>] = J dx-T[4>] , (211) 
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is of the form (209), because the integration can be viewed as a projection onto 
the longitudinal momentum sector k + = 0. The most important examples for 
such quantities are the Poincare generators, as is obvious from the representa- 
tions (107, 108). This implies in particular that the trivial light-cone vacuum is 
an eigenstate of the fully interacting light-cone Hamiltonian P~ , 

P~\0) = (0|P-|0}|0) . (212) 

This can be seen alternatively by considering 

P+P-|0) =P-P + \0) =0, (213) 

which says that P~\0) is a state with fc+ = 0, so that P~ must have a Fock 
representation like U in (209). 

The actual value of (0|P~|0) is not important at this point as it only defines 
the zero of light-cone energy. Note that, within the instant form, the Fock or 
trivial vacuum is not an eigenstate of the full Hamiltonian as the latter usually 
contains terms with only creation operators where positive and negative three- 
momenta compensate to zero as in (207). The instant-form vacuum thus is 
unstable under time evolution. Such a vacuum, a typical example of which is 
provided by (207), is called 'nontrivial'. 

This concludes the general discussion of light-cone quantum field theory. 
Having understood the foundations of the approach we are now heading for the 
applications. We begin with a survey of the light-cone Schrodinger equation. 
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4 Light-Cone Wave Functions 

In this section we collect some basic facts about the eigenvalue problem of the 
light-cone Hamiltonian, or, in other words, about the light-cone Schrodingcr 
equation and its solutions, the light-cone wave functions. Throughout this sec- 
tion, I will use the conventions of (Brodsky and Lepage 1989). 

4.1 Kinematics 

To set the stage for the definition of light cone wave functions let me first intro- 
duce some relevant kinematical variables. Consider a system of many particles 
which, for the time being, will be assumed as noninteracting. Let the i th particle 
have mass m, and light-cone four-momentum 

Pi = (pt,P±i,Pi) ■ (214) 

As the particles are free, the total four-momentum is conserved and thus given 
by the sum of the individual momenta, 

P = J>i- (215) 

i 

The individual particles are on-shell, so their four-momentum squared is 

Pi=PiP;-p 2 Li = ™ 2 i ■ (216) 

The square of the total four-momentum, on the other hand, defines the free 
invariant mass squared, 

P 2 = p+p - _pl= Ml , (217) 

a quantity that will become important later on. We introduce relative momentum 
coordinates Xi and k±i via 

Pt = x t P + , (218) 
p ±4 = £ 4 P ± + k ±t . (219) 

Thus, Xi and kj_, denote the longitudinal momentum fraction and the relative 
transverse momentum of the i th particle, respectively. Comparing with (215) we 
note that these variables have to obey the constraints 

5>j = l, ^k ±4 = 0. (220) 

i i 

A particularly important property of the relative momenta is their boost invari- 
ance. To show this we calculate, using (86), 



x'i = e"pt/e"P + = Xi 



(221) 
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From this and (88) we find in addition 

k± 4 = P±i - ^ P '± = P± 4 + - ^( p ^ + v ^ p+ ) = k ^ - ( 222 ) 

which indeed proves the frame independence of Xi and k±i. 

Let us calculate the total light-cone energy of the system in terms of the 
relative coordinates. Making use of the constraints (220), we obtain 

— Zr-^ . — I — m? \ 

P CM +P-. (223) 

This is another important result: the light-cone Hamiltonian P~ separates into 
a center-of-mass term, 

P6m=Pl/P + , (224) 
and a term containing only the relative coordinates, 

(225, 

The second identity, which states that P~ is essentially the free invariant mass 
squared, follows upon multiplying (223) by P+, 

P + P- = P + P- - P| = M 2 = + m ' ■ (226) 

% 

To simplify things even more, one often goes to the 'transverse rest frame' where 
and therefore the center-of-mass Hamiltonian Pq M from (224) vanish. 
In the interacting case, the dynamical Poincare generators acquire 'potential' 
terms as I have shown in Section 3.1. The light-cone Hamiltonian, e.g. becomes 
P~ = P ~ + V, leading to a four-momentum squared 

P 2 = P + (P Q +V) - Pl = M 2 . (227) 

Subtracting (217) we obtain the useful relation, 

M 2 — Mq = P + V = W . (228) 

In the quantum theory, this operator identity, when applied to physical states, 
is nothing but the light-cone Schrodinger equation. 

Summarizing we note that the special behavior under boosts together with 
the transverse Galilei invariance leads to frame independent relative coordinates 
and a separation of the center-of-mass motion, reminiscent of ordinary nonrelati- 
vistic physics. This is at variance with the instant form, where the appearance of 
the notorious square root in the energy, P° = (P 2 + M^) 1 / 2 , prohibits a similar 
separation of variables. 
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4.2 Definition of Light Cone Wave Functions 

Let us first stick to a discrete notation and, for the time being, stay in 1+1 
dimensions. We thus have a Fock basis of states 

\n) = 4|0) , 
\m,n) = 4,410) , 



|n 1 ,...,n iV } = 4 1 ---4 N |0} ■ (229) 
This leads to a completeness relation defining the unit operator in Fock space, 

1 = |0)(0| + ^|n)(n| + i E \m,n)(m,n\ + ... 

n>0 m,n>0 
1 

= |0}<0|+ E ^ E \ni,...,n N )(n u ...,n N \. (230) 

JV=1 ' ni,...,njv>0 

An arbitrary state can thus be expanded as 

|V> = £W>l«> + 5 E (™,n\TP)\m,n) + ... . (231) 

n>0 m,n>0 

The sums are such that the longitudinal momenta in each Fock sector add up to 
the total longitudinal momentum of \tp). Note that the vacuum state does not 
contribute as it is orthogonal to any particle state, (Q\ip) = 0. The normalization 
of this state is obtained as 

(V#)=]T|(n|^| 2 + ± Yl \(^n\^ + ... 

n>0 m.n>0 

oo 

= Ei E \(m,...,n N m 2 . (232) 

N=l ' ni,...,n N >0 

Let us assume that the state \ip) corresponds to a bound state obeying the 
light-cone Schrodinger equation derived from (228), 

(M 2 — Mq)\iP) — W\tp) . (233) 

We want to project this equation onto the different Fock sectors. For this we need 
the eigenvalues of the free invariant mass squared when applied to an A^-particle 
state 

\N) = \n u ...n N ) . (234) 
We find that \N) is an eigenstate of M$ , 

N 2 

M%\N)=Y—\N)^M%\N). (235) 

i—1 
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The light-cone Schrddinger equation thus becomes a system of coupled eigen- 
value equations, 



(m 2 - Mi) m 

(M 2 - Mf) (klty) 



(l\W\m) (l\W\mn) 
(kl\W\m) (kl\W\mn) 



(236) 



Clearly, this represents an infinite number of equations which in general will 
prove impossible to solve unless the matrix is very sparse and/or the matrix 
elements are small. The former condition is usually fulfilled as the interaction 
W in the light-cone Hamiltonian typically changes particle number at most 
by two 7 . Assuming the matrix elements to be small amounts to dealing with 
a perturbative situation. This will be true for nonrelativistic bound states of 
heavy constituents, but not for light hadrons which we are mainly interested in. 
There are, however, situations where the magnitude of the amplitudes decreases 
enormously with the particle number TV, so that it is a good approximation to 
restrict to the lowest Fock sectors. In instant form field theory this has long been 
known as the Tamm-Dancoff method (Tamm 1945; Dancoff 1950). 

Let us turn to the more realistic case of 3+1 dimensions in a continuum 
formulation. The invariant normalization of a momentum eigenstate \P + , Pj_) = 
\P) is given by 

(P\K) = 16tt 3 P+5 3 (P - K) . (237) 

We already know that the bare Fock vacuum is an eigenstate of the interacting 
Hamiltonian. It thus serves as an appropriate ground state on top of which we 
can build a reasonable Fock expansion. If we specialize immediately to the case 
of QCD, we are left with the Fock basis states 



10} , 

\qq : fej.aj) = 6 f (fei, ai)d t (fc 2 , a 2 )|0) , 

\qqg : k i} ai) = 6 f (fci, a 1 )d^(k 2 , a 2 )a t (fc 3 , a 3 )|0) , 



(238) 



(239) 



In these expressions, b\ d) and a) create quarks q 1 antiquarks q and gluons g 
with momenta ki from the trivial vacuum |0). The a, denote all other relevant 
quantum numbers, like helicity, polarization, flavor and color. 

In a more condensed notation we can thus describe, say, a pion with momen- 
tum P = (P+,P_l), as 

JJdxi 16 J 3 V n/7r(^ i k_Li, Xi) n : XiP + , #iP_i_+kxi, Xi^ , (240) 



7 Note that terms with only creation operators are forbidden by fc + -conservation. 
Still, in 1+1 dimensions, things can become messy as interactions with polynomials 
of arbitrary powers in 4> are allowed without spoiling renormalizability (Zinn- Justin 
1996) 
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where we have suppressed all discrete quantum numbers apart from the helicities 
A, . The integration measure takes care of the constraints (220) which the relative 
momenta in each Fock state (labeled by n) have to obey, 



(241) 
(242) 



As a mnemonic rule, we note that any measure factor d 2 k±i in (240) is always 
accompanied by 1/167T 3 . 

The most important quantities in (240) are the light-cone wave functions 



Y[dxi = \\ dx l s ( i - 

i i j 

\\d 2 ku = 16^ 3 I] d 2 k ±l S 2 ( k ±J 



^ n / 7: (x l ,k_ Ll ,X i ) = (n : x l P + ,x l P ± +\s.± i ,X i \ir(P)) , 



(243) 



which are the amplitudes to find n constituents with relative momenta fi\ = 
XiP + , p ±i — XiP± + k±i and helicities Ai in the pion. Due to the separation 
properties of the light-cone Hamiltonian the wave functions do not depend on 
the total momentum P of the pion. Applying (237) to the pion state (240), we 
obtain the normalization condition 



(244) 



71, Xi 



The light-cone bound-state equation for the pion is a straightforward general- 
ization of (236), 



(M 2 - M 2 ) (qq\Tr) 
(M 2 - M% g ) (qqg\n) 



(qq\W\qq) (qq\W\qqg) . . . 
(qqg\W\qq) (qqg\W\qqg) . . . 



[qq\v 
(qqgW) 



(245) 



If a constituent picture for the pion were true, the valence state would dominate, 

IV>2A| 2 » KaI 2 . ™> 2 > ( 246 ) 

and, in the extreme case, the pion wave function would be entirely given by the 
projection (qq\n) onto the valence state. All the higher Fock contributions would 
vanish and the unitarity sum (244) would simply reduce to 

Y.[ dx j ^l^(^AA')| 2 = l. (247) 

A A' 

We will later discuss a model where this is indeed a good approximation to 
reality. 
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4.3 Properties of Light Cone Wave Functions 

Let us rewrite the light-cone bound-state equation (236) by collecting all light- 
cone wave functions ip n = (n\ip) into a vector 

* = m=ns- (248) 

From this expression it is obvious that all light-cone wave functions tend to 
vanish whenever the denominator 

e = M 2 - M 2 = M 2 - (Y,PiY = m2 - E k2±i+m ' (249) 

becomes very large. This quantity measures how far off energy shell the total 
system, i.e. the bound state is, 

p--£>r= e /P+. (250) 

i 

For this reason, e is sometimes called the 'off-shellness' (Namyslowski 1985; 
Lavelle et al. 1987). We thus learn from (248) that there is only a small overlap 
of the bound state with Fock states that are far off shell. This implies the limiting 
behavior 

V>(zj,k_u, Xi) — > for Xi — > , k\ { — > oo . (251) 

These boundary conditions are related to the self-adjointness of the light-cone 
Hamiltonian and to the finitcness of its matrix elements. Analogous criteria have 
been used recently to relate wave functions of different Fock states n (Antonuccio 
et al. 1997) and to analyse the divergence structure of light-cone perturbation 
theory (Burkardt 1998). 

Omitting spin, flavor and color degrees of freedom, a light-cone wave function 
will be a scalar function 0(xj, kj_j) of the parameter e. This is used for building 
models, the most common one being to assume a Gaussian behavior, originally 
suggested by (Terent'ev 1976), 

0(x 4 ,k ±4 )=iVcxp(-|e|//3 2 ), (252) 

where f3 measures the size of the wave function in momentum space. Note, how- 
ever, that a Gaussian ansatz is in conflict with perturbation theory which is the 
appropriate tool to study the high behavior and indicates a power decay 
of the renormalized wave function — up to possible logarithms (Brodsky and 
Lepage 1989). For the unrenormalized wave functions the boundary conditions 
(251) are violated unless one uses a cutoff as a regulator (see next section). 

As the off-shellness e is the most important quantity characterizing a light- 
cone wave function let us have a closer look by specializing to the simplest 
possible system, namely two bound particles of equal mass. One can think of 
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this, for instance, as the valence wave function of the pion or positronium. The 
off-shellness becomes 

_ M 2 k\+m 2 _ 1 r , ,o / 1x2 4m 2 - M 2 



x(l — x) x(l — x) 



M 2 {x-\Y + +fci . (253) 



>o 



The second term in square brackets is positive because, for a bound state, the 
binding energy, 

E = M — 2m , (254) 

is negative so that 2m > M. As a result, the off-shellness is always negative. 
Only for free particles it is zero, because all momentum components (including 
the energy) sum up to the total momentum. In this case, each individual term 
in (253) vanishes, 

M = 2m , x=\, kj_ = . (255) 

It follows that the light-cone wave function of a two-particle system (composed 
of equal-mass constituents) with the binding energy switched off 'adiabatically', 
is of the form 

<j)(x,k ± ) ~ 5(x- l/2)<5 2 (k ± ) . (256) 



4.4 Examples of Light Cone Wave Functions 

From the discussion above, one expects that for weak binding, in particular for 
nonrelativistic systems, the wave functions will be highly peaked around x = 1/2 
(in the equal mass case) and kj^ = 0. Let us check this explicitly for hydrogen 
like systems which constitute our first example (Lepage et al. 1981). 



Example 1: Hydrogen— like Systems. Let me recall the ordinary Schrodinger 
equation of the Coulomb problem written in momentum space, 

{ E - Q «*> - / - - / ■ <-> 

This integral equation looks very similar to a light-cone Schrodinger equation 
within a two-particle truncation. The Coulomb kernel is due to the exchange 
of an instantaneous photon having a propagator proportional to 5{x°). One can 
actually solve the Coulomb problem directly in momentum space (Fock 1935; 
Bethe and Salpeter 1957) but for our purposes it is simpler just to Fourier 
transform the ground state wave function VoM = N exp(— mar), yielding 

ma 

MP) = SttJV , (258) 

(p^ + m z a z ) z 

with a = e 2 /Air = 1/137 being the fine structure constant, m the reduced mass 
and p the relative momentum. 
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How does this translate into the light-cone language? To answer this ques- 
tion, we go to the particle rest frame with P = or P + = P~ = M and 
P_i_ = 0, implying p ±i = kj_j. In this frame, the nonrelativistic limit is de- 
fined by the following inequalities for the constituent masses and momenta (in 
ordinary instant-form coordinates), 

^-m l ~^«|p 4 |«m l . (259) 

The prototype systems in this class are of hydrogen type where we have for 
binding energy and r.m.s. momentum, 

i*i - g - . m 

(p) = ma . (261) 
In this case, the hierarchy (259) becomes 

a 2 

— < a < 1 , (262) 

which is fulfilled to a very good extent in view of the smallness of a. 
Consider now the longitudinal momentum of the i th constituent, 

p+ = p° + p» ~ mi + £j- + p\ = Xi P + = Xl M . (263) 

We thus find that we should replace p\ in instant-form nonrelativistic wave 
functions by 

p\ = x t M - m, , (264) 

where we neglect terms of order p?/mj. Let us analyze the consequences for the 
off-shellness. The latter is in ordinary coordinates 

e = M 2 -M 2 = (M + J>°)(M-£>°). (265) 

i i 

We thus need 

Erf-E- + E^- = «-^E^-. <*»> 

tit i 

with E — M — J2 m i denoting the mass-defect, which is a small quantity, E <C 
M. The off-shellness (265), therefore, becomes 

e * 2M ( E - £ p? /2ro< ) * -2M £ ^ + ^ " ^ , (267) 
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where we have performed the replacement (264) in the second identity. The 
light cone wave functions will be peaked where the off-shcllness is small, that 
is, for 

Xi =mi/M, and k_u = , (268) 

as expected from the noninteracting case. 

To be explicit, we consider the ground state wave function of positronium, 
given by (258) with the reduced mass m being half the electron mass m e . Using 
the replacement prescription (264) once more, we obtain 

171 d 

4>(x, k ± ) = 8irN ^ , (269) 

k\ + (xM — m e ) 2 + m 2 a 2 

where M ~ 2m e is the bound state mass. This result is valid for small momenta, 
i.e. when k\, (xM — m e ) 2 -C m 2 . It is obvious from (269) that the positronium 
wave function is sharply peaked around x = m e /M ~ 1/2 and k\ = 0. 



Example 2: 't Hooft Model. The 't Hooft model ('t Hooft 1974; 't Hooft 
1975) is QCD in two space-time dimensions with the number Nc of colors being 
infinite. The Lagrangian is 

£ = - m )i> - \ F ^ F>XV • (270) 

The limit of large Nc is taken in such a way that the expression g 2 Nc, g denoting 
the Yang-Mills coupling, stays finite. In two dimensions, g has mass dimension 
one, which renders the theory superrenormalizable and provides a basic unit of 
mass, namely, 

Mo = ^g 2 N c /2ir . (271) 

The model is interesting because it contains physics analogous or similar to what 
one finds in 'real' QCD 8 . First of all, the model is (trivially) confining due to the 
linear rise of the Coulomb potential in 2d. This is most easily exhibited by work- 
ing in light-cone gauge, A + = 0, and eliminating A~ via Gauss's law. In this way 
it becomes manifest that there are no dynamical gluons in 2d. Within covariant 
perturbation theory, the Coulomb potential can be understood as resulting from 
the exchange of an instantaneous gluon. 

In the next section, we will discuss the spontaneous breakdown of chiral sym- 
metry in QCD. It turns out that in the 't Hooft model a similar phenomenon 
occurs: chiral symmetry is 'almost' spontaneously broken (Witten 1978; Zhitnit- 
sky 1986). As a consequence, there arises a massless bound state in the chiral 
limit of vanishing quark mass ('t Hooft 1974; 't Hooft 1975) which we will call 



For a recent review on the 't Hooft model, see (Abdalla and Abdalla 1996). 
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the 'pion' for brevity 9 . Furthermore, there is a nonvanishing quark condensate 
in the model which has first been calculated by (Zhitnitsky 1986), 

(0|V^|0) /N c = -0.28868 . (272) 

Note that the condensate is proportional to Nq as it involves a color trace. 

It turns out that the 't Hooft model has one big advantage for the application 
of light-cone techniques which is due to the large-iVc limit. The matrix elements 
entering the light-cone Schrodingcr equation in the two-particle sector have the 
following Afc-dependence, 

( W n)~(£j^) ~N c n . (273) 

As a result, those diagrams which correspond to a change in particle number 
(like 2 — ► 4, 2 — > 6, . . .) arc suppressed by powers of 1/Nq- The truncation to the 
two-particle sector therefore becomes exact: the 'pion' is a pure quark-antiquark 
state; there are no admixtures of higher Fock states. A constituent picture is thus 
realized, and we are left with the light-cone Schrodingcr equation, 



M° m 



x(l — x) 



This expression defines the 'Coulomb problem' of the 't Hooft model. It corre- 
sponds to the first line of (245) where, as a result of the truncation, only the 
matrix element W% = (qq\W\qq) has been retained. We will refer to (274) as the 
't Hooft equation in what follows. <j>(x) denotes the valence part of the 'pion' 
wave function, x and y are the momentum fractions of the two quarks (with 
equal mass m) in the meson. The symbol V indicates that the integral is defined 
as a principal value ft Hooft 1974; Gelfand and Shilov 1964). It regularizes 10 
the Coulomb singularity l/(x — y) 2 in the matrix element W^- All masses are 
expressed as multiples of the basic scale fio defined in (271). The eigenvalue M 
denotes the mass of the lowest lying bound state (the 'pion'). Our objective is 
to calculate M and <p. 

In his original work on the subject, ('t Hooft 1974; 't Hooft 1975), 't Hooft 
used the following ansatz for the wave function 

<j){x) = x^{l-xf. (275) 

This ansatz is symmetric in x <-> 1 — x (charge conjugation odd), and (i is sup- 
posed to lie between zero and one so that the endpoint behavior is nonanalytic. 



As explained in (Witten 1978; Zhitnitsky 1986), this is not in contradiction with 
Coleman's theorem (Coleman 1973) as the 'pion' is not a Goldstone boson. 
(Wu 1977) has suggested a theoretical alternative to the principal value which nowa- 
days is called 'Leibbrandt-Mandelstam prescription (Mandelstam 1983; Leibbrandt 
1984). It leads to completely different physics. This apparent contradiction has only 
recently been clarified (Bassetto et al. 2000). 
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As a nontrivial boundary condition, one has the exact solution of the masslcss 
case, to = (the 'chiral limit'), 

M 2 = 0, and <j)(x) = l, i.e. (3 = 0. (276) 

In this limit, the 'pion' wave function is constant, i.e. the pion has no structure 
and is point-like. The main effect of a nonzero quark mass is the vanishing of 
the wave functions at the endpoints implying a nonzero (3. This suggests the 
following series expansion for (3, 

(3{m) = (3im + f3 2 m 2 + (3 s m 3 + (9(m 4 ) , (277) 

and for the 'pion' mass squared, 

M 2 = Mim + M 2 to 2 + M 3 m 3 + 0(to 4 ) . (278) 

As is obvious from the last two expressions, we are working in the limit of small 
quark mass to. The expansion (277) shows that also (3 is small in this case so 
that the wave function will be rather flat (for intermediate values of x). On the 
other hand we know from Example 1 that light-cone wave functions are highly 
peaked near x = 1/2 in case the binding is weak. This suggests that, for small 
quark mass to, the 'pion' is rather strongly bound. 

The exponent f3 in (275) can actually be determined exactly by studying the 
small- 2 behavior of the bound state equation (274). To this end we evaluate 
the principal value integral for x — > and plug it into (274). This yields the 
transcendental equation ft Hooft 1974), 

m 2 - 1 + 7T/3 cot 7T/3 = . (279) 

Using this expression we can determine [3 either numerically for arbitrary m or 
analytically for small m, which yields the coefficients of (277), 

0=— to + 0(to 3 ). (280) 

The 'pion' mass is determined by calculating the expectation value of the light- 
cone Hamiltonian (274) in the state given by 't Hooft's ansatz (275) 11 . This 
yields M 2 as a function of (3 and to, 

M 2 = J-to + 0(to 2 ) . (281) 
Pi 

Upon comparing with (277) and (280) the lowest order coefficient Mi, which is 
the slope of M 2 (m) at to = 0, is found to be 

Mi = 2tt/V3 . (282) 

Note that M 2 indeed vanishes in the chiral limit. 



The relevant integrals can be found in (Bardeen et al. 1980) and (Harada et al. 1994). 
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Having obtained an approximate solution for the mass and wave function of 
the 'pion' we are in the position to calculate 'observables'. It turns out that, for 
small m, all of them can be expressed in terms of the lowest order coefficient, 
Mi. We begin with the 'pion decay constant' (Callan et al. 1976; Zhitnitsky 
1986), which is given by the 'wave function at the origin', i.e. the integral over 
the (momentum space) wave function 

iNr M 2 f iNr 

U = (o|v>H 5 V>k) = \ — 7T- / dx ^ x ) = \hr M i • ( 283 ) 

V 7r Ira J V 47r 

The quark condensate is obtained via a sum rule using the chiral Ward identity 
(Zhitnitsky 1986; Heinzl 1996b), 

(0|^|0> = -mjfc = -_£ Ml . (284) 

Inserting the value (282) for Mi this coincides with (272). The last identity can 
actually be viewed as the 'Gell-Mann-Oakes-Renner relation' (Gell-Mann et al. 
1968) of the 't Hooft model, 

M 2 = — 47r ra (0|^V|0) /Nc + 0(m 2 ) , (285) 

which provides a relation between the particle spectrum (the 'pion' mass) and a 
ground state property (the condensate). This is conceptually important because 
it implies that we can circumvent the explicit construction of a nontrivial vacuum 
state by calculating the spectrum of excited states, i.e. by solving the light- 
cone Schrodinger equation. The eigenvalues and wave functions actually contain 
information about the structure of the vacuum! This point of view has been 
adopted long ago in the context of chiral symmetry breaking within the (light- 
cone) parton model (Casher et al. 1971): "In this framework the spontaneous 
symmetry breakdown must be attributed to the properties of the hadron's wave 
function and not to the vacuum" (Casher and Susskind 1974). Related ideas 
have been put forward more recently in (Lenz et al. 1991). 

In the above, we have been using the value for (3 given in (280). One can 
equally well use /3 as a variational parameter and minimize the expectation 
value of the mass operator with respect to it. This yields the same result for Mi, 
namely (282). The variational method, however, is better suited if one wants to 
go beyond the leading order in expansion (278). This has been done in (Harada 
et al. 1998). The results are shown in Table 2 where we list the expansion coeffi- 
cients Mi of the 'pion' mass squared, M 2 , as they change with increasing number 
of variational parameters, (a, b, c, d). 

Interestingly, the value of Mi does not change at all by enlarging the space of 
trial functions. M 2 and M 3 , on the other hand, do change and show rather good 
convergence. For M 2 we finally have a relative accuracy of 8- 1CT 7 , and for M 3 of 
4 • 10~ 5 . Furthermore, the coefficients are getting smaller if one adds more basis 
functions, in accordance with the variational principle. The associated light-cone 
wave functions are shown in Fig. 8. There are only minor changes upon including 
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more variational parameters. In Fig. 9 we display the eight lowest excited states. 
They habe been obtained using the position of the nodes as additional variational 
parameters (Stern 1999). 



Table 2. Expansion coefficients of M 2 for the 't Hooft model obtained by successively 
enlarging the space of variational parameters. Mi is the standard 't Hooft result (282). 
Note the good convergence towards the bottom of the table. 



Ansatz Mi = 2ir/V3 M 2 M 3 

't Hooft 3.62759873 3.61542218 0.043597197 

a 3.62759873 3.58136872 0.061736701 

b 3.62759873 3.58107780 0.061805257 

c 3.62759873 3.58105821 0.061795547 

d 3.62759873 3.58105532 0.061793082 



Fig. 8. The light-cone wave function of the 't Hooft model 'pion' for m = 0.1. The 
solid curve represents the result from 't Hooft's original ansatz, the dashed curve our 
best result (with maximum number of variational parameters). At the given resolution, 
however, the curves of all extensions of 't Hooft's ansatz (a,b,c,d) lie on top of each 
other. 
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Fig. 9. Wave functions of the first eight excited states in the 't Hooft model as obtained 
via variational methods (Stern 1999). Note that all wave functions vanish at the end 
points, x = and x = 1. 
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Example 3: Gaussian Model. A very simple and intuitive example of a light- 
cone wave function is provided by the Gaussian model in the form presented by 
(Lepage et al. 1981), 



(j)(x,k±) = N exp 



1 k\ + m 2 



A 2 x(l - x) 



(286) 



This is a scalar two-particle wave function that drops exponentially with the 
free invariant mass squared of two constituents with equal mass m. The short- 
coming of the model is that the wave function is not derived dynamically i.e. as 
a solution of a light-cone Schrodinger equation. However, it satisfies all the re- 
quirements derived in Subsection 4.3. The benefit of the model is its simplicity. 
The latter is even enhanced if one neglects the constituent mass m by assuming 
m <C A. Then one is left with only two parameters, the normalization N and 
the transverse size A. 

One enforces a constituent picture by normalizing to unity, 

f 1 f H 2 k, N 2 A 2 i 

U(x,^)\\ 2 ^J Q dxJ ^\x M = — 2 = 1 . (287) 

This obviously relates N and A. If we now calculate the r.m.s. transverse mo- 
mentum in the bound state described by (286), we find using (287), 



' ' / / J^ k '* 2 ^ ) = 4^ = A ' /m ' (288) 



If we view the Gaussian wave function as a crude model for, say, the pion we 
can actually estimate the width parameter A. The pion is highly rclativistic, so 
we expect its r.m.s. transverse momentum to be of the order of the constituent 
quark mass, (fc^) 1 / 2 — m,Q ~ 330 MeV. This leads to a typical width of A ~ 1 
GeV. 

Having thus determined the two parameters of the model we could go on and 
calculate observables (Lepage et al. 1981). This will actually be done in the next 
section using a more realistic model for the (determination of) the pion wave 
function. 



5 The Pion Wave Function in the NJL Model 

The ultimate goal of light-cone field quantization is to derive and solve light- 
cone Schrodinger equations, in particular the one of QCD. This would yield 
hadron masses and light-cone wave functions and thus detailed information on 
the internal structure of mesons and baryons. In order to successfully pursue 
this program, a number of problems has to be overcome. 

In dealing with gauge theories in a Hamiltonian framework one has to solve 
Gauss's law together with the light-cone specific constraints. The only known 
solution so far is in the light-cone gauge A + = (Tomboulis 1973; Casher 1976; 



60 Thomas Heinzl 



Lepage et al. 1981) which, however, is beset by infrared problems of its own — see 
e.g. (Bassetto et al. 1991) 12 . A Hamiltonian formulation analogous to the Weyl 
gauge, where one fixes the residual gauge freedom and solves Gauss's law after 
quantization, seems particularly difficult (Heinzl 1996a). 

Even after successful light-cone quantization of QCD one encounters a se- 
vere problem: the theory has to be renormalized. Otherwise, light-cone wave 
function will not be normalizable and physically meaningful. Rcnormalization in 
a Hamiltonian framework presents enormous difficulties as there is no explicit 
covariancc. In addition, rotational invariance is not manifest within the front 
form (see Section 2). Due to the lack of these important symmetries, there is an 
abundance of possible counterterms which even can be nonlocal, e.g. behave like 
<~ l/fc + . As a result, to the best of my knowledge, the renormalization program 
has not been extended beyond one loop — with one notable exception in QED 
(Brodsky et al. 1973). 

The issue of renormalization is, of course, delicately intertwined with solv- 
ing the light-cone bound state equation. The latter attempt will in general 
only be feasible if a Tamm-Dancoff truncation in particle number is performed 
(Harindranath et al. 1990). This again violates important symmetries (even 
gauge symmetry). One hopes, however, that a Wilsonian renormalization group 
explicitly taylored for this case will restore these symmetries (Wilson et al. 1994). 
The present status of this program is nicely reviewed in (Perry 1999). 

A conceptual problem has already been mentioned. The light-cone vacuum is 
trivial as argued in Section 3. On the other hand, it is well known that the instant 
form vacuum is populated by all sorts of nonperturbative quantum fluctuations 
leading to nonvanishing vacuum expectation values or condensates. The whole 
business of QCD sum rules (Shifman et al. 1979) is based upon this picture. How 
do we reconcile this with the triviality of the light-cone vacuum? A possible 
resolution to this 'triviality problem' has been given in the last section. The 
spectrum of excited states actually carries implicit information on the structure 
of the vacuum. The task then is to make this information explicit. 

The list of problems just given is yet another manifestation of the 'principle 
of conservation of difficulties'. In a first attempt to tackle these problems I will 
simply side-step most of them by considering an (instructive) model instead 
of QCD. This model, however, is designed to capture some important physical 
features of 'real' QCD. The idea is originally due to Nambu and Jona-Lasinio 
(NJL), who, back in 1961, invented a "Dynamical Model of Elementary Particles 
Based on an Analogy with Superconductivity" (Nambu and Jona-Lasinio 1961). 
It was meant to provide a microscopic mechanism for the generation of nucleon 
masses, with the mass gap being the analog of the BCS energy gap in a super- 
conductor. Nowadays, with the nucleons replaced by quarks, the model serves as 
a low-energy effective theory of QCD explaining the spontaneous breakdown of 
chiral symmetry. Let me thus give a brief introduction to the latter phenomenon 



Recently, some progress has also been made in the covariant gauge (Srivastava and 
Brodsky 2000). 
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before I come to the detailed explanation of the model. 



5.1 A Primer on Spontaneous Chiral Symmetry Breaking 

If we have a look at Table 3, which provides a list of all quark flavors, we realize 
that there are large differences in the quark masses as they appear in the QCD 
Hamiltonian. In particular, there is a hierarchy, 

m u ,m d <m s < m c ,mb,m t . (289) 

V v ' V v ' 

light quarks heavy quarks 

As the masses of heavy and light quarks are separated by the very same scale 
(~ 1 GcV) as the perturbative and nonperturbative regime, one expects differ- 
ent physics associated with those two kinds of quarks. This expectation turns 
out to be true. The physics of heavy quarks is governed by a symmetry called 
'heavy quark symmetry' leading to a very successful 'heavy quark effective the- 
ory' (Neubert 1994). The physics of light quarks, on the other hand, is governed 
by chiral symmetry which we are now going to explain. 



Table 3. The presently observed quark Savors. Q/e is the electric charge in units of 
the electron charge. The (scale dependent!) quark masses are given for a scale of 1 
GeV. 



flavor 




Q/e 


mass 




down 


d 


-1/3 


10 


MeV 


up 


u 


+2/3 


5 


MeV 


strange 


s 


-1/3 


150 


MeV 


charm 


c 


+2/3 


1.5 


GeV 


bottom 


b 


-1/3 


5.1 


GeV 


top 


t 


+2/3 


180 


GeV 



Let us write the QCD Hamiltonian in the following way, 

Hqcb = H x + ^M^j , (290) 

A4 = diag(m u ,md,m s ) being the mass matrix for the light flavors. To a good 
approximation, one can set M. = 0. In this case, the QCD Hamiltonian H x is 
invariant under the symmetry group SU(3)r <g> SU(3)l, the chiral flavor group. 
Under the action of this group, the left and right handed quarks independently 
undergo a chiral rotation. Due to Noether's theorem, there are sixteen conserved 
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quantities, eight vector charges and, more important for us, eight pscudo-scalars, 
the chiral charges QS? satisfying 

[QS,H x ]=0. (291) 

This states both that the chiral charges are conserved, and that H x is chirally 
invariant. Under parity, — ► —Q%. Now, if \A) is an eigenstate of H x , so is 
Qgl-A) with the same eigenvalue. Thus, one expects (nearly) degenerate parity 
doublets in nature, which, however, do not exist empirically. The only explana- 
tion for this phenomenon is that chiral symmetry is spontaneously broken. In 
contradistinction to the Hamiltonian, the QCD ground state (the vacuum) is not 
chirally invariant, 

Qg|0) ^ . (292) 

For this reason, there must exist a nonvanishing vacuum expectation value, the 
quark condensate, 

(^V) = {Ml + Mr) • (293) 

This condensate is not invariant (it mixes left and right) and therefore serves 
as an order parameter of the symmetry breaking. Note that in QCD the quark 
condensate is a renormalization scale dependent quantity. A recent estimate can 
be found in (Dosch and Narison 1998), with a numerical value, 

(0|^V|0)(1 GeV) ~ (-229 MeV) 3 . (294) 

The spontaneous breakdown of chiral symmetry thus implies that the (QCD) 
vacuum is nontrivial: it must contain quark-antiquark pairs with spins and mo- 
menta aligned in a way consistent with vacuum quantum numbers. A possible 
analog is the BCS ground state given in (207). 
In terms of the full quark propagator, 

^ - ^m> • 

where we have allowed for a momentum dependent (or 'running') mass M(p), 
the quark condensate is given by 

*.,/^.*/^. (2M> 

We thus see that the involved Dirac trace yields a nonvanishing condensate 
only if the effective quark mass M(p) is nonzero. This links the existence of 
a quark condensate to the mechanism of dynamical mass generation. In this 
way we have found another argument that the bare quarks appearing in the 
QCD Hamiltonian (290) indeed acquire constituent masses. Of course we are 
still lacking a microscopic mechanism for that. 

Goldstone's theorem (Goldstone 1961) now states that for any symmetry 
generator which does not leave the vacuum invariant, there must exist a massless 
boson with the quantum numbers of this generator. This results in the prediction 
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Table 4. Masses of the pseudoscalar octet mesons (in MeV). 



meson 






K°, K° 




V 


mass 


135 


140 


500 


494 


549 



that in masslcss QCD one should have an octet of massless pseudoscalar mesons. 
In reality one finds what is listed in Table 4. 

The nonvanishing masses of these mesons are interpreted as stemming from 
the nonvanishing quark masses in the QCD Hamiltonian which break chiral 
symmetry explicitly Being small, they can be treated as perturbations. 

It should be stressed that chiral symmetry has nothing to say about the 
mechanism of confinement which presumably is a totally different story. This 
is also reflected within a QCD based derivation of chiral symmetry breaking in 
terms of the instanton model, as presented e.g. in (Diakonov 1996). This model 
explains many facts of low-energy hadronic physics but is known not to yield 
confinement. It is therefore possible that confinement is not particularly relevant 
for the understanding of hadron structure (Diakonov 1996). 

The instanton vacuum actually leads to an effective theory very close to the 
NJL model, the basics of which arc our next topic. 

5.2 NJL Folklore 

Before I consider the light-cone formulation of the model, let me briefly recall 
its main physical features 13 . In its standard form, the NJL model has a chirally 
invariant four-fermion interaction, which can be imagined as the result of 'in- 
tegrating out' the gluons in the QCD Lagrangian. For simplicity I concentrate 
on the case of one flavor. Extension to several flavors is straightforward. In the 
chiral limit (quark mass mo = 0), the Lagrangian is 

C = $%frl> - G^^f = Co + £ int . (297) 

Its four-fermion interaction is chirally symmetric under U(1)l x U (1)_r- Wc shall 
see in a moment that this symmetry is spontaneously broken. 

It is important to observe that the coupling G of the model has negative mass 
dimension, [G] = —2, hence, it is not renormalizable. Accordingly, it requires a 
cutoff which is viewed as a parameter of the model that is to be fixed by phe- 
nomenology. We thus follow the general spirit of effective field theory (Holstcin 
1995; Kaplan 1995; Manohar 1996). From the point of view of the light-cone 



For recent reviews, see (Vogl and Weise 1991; Klevansky 1992; Hatsuda and Kunihiro 
1994). 
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formulation to be developed later, the nonrenormalizability is an advantage: it 
enables us to 'circumvent' the difficulties of the light-cone renormalization pro- 
gram. There simply is no 'need' to renormalize. 

Following the standard approach (Nambu and Jona-Lasinio 1961; Vogl and 
Weise 1991) wc treat the model in mean-field approximation (which actually 
coincides with the large-iV limit). We begin with a Fierz transformation by 
schematically rewriting the interaction Lagrangian in Fierz symmetric form 

Ant = G^c^V) 2 - G(^) 2 + . . . , (298) 

i 

where i = S,P,V, A enumerates the different Dirac bilinears (scalar, pseu- 
doscalar, vector and axial vector, respectively). In (298), I have only displayed 
the scalar part of the interaction because only the scalar density S = vpip can 
have a vacuum expectation value, the quark condensate, 

(S) = (0|^|0) . (299) 

Let us determine this quantity in mean-field approximation. To define the latter 
we calculate, 

S 2 = (S - (S) + (S)) 2 = (S- (S)) 2 + 2S(S) - (S) 2 ~ 2S(S) + const . (300) 

Thus, by neglecting quadratic fluctuations of S around its expectation value, we 
linearize the interaction and obtain the mean-field Lagrangian, 

£mfa = ^(^ + 2G(5))V. (301) 

The mean-field solution has a very intuitive explanation. One essentially argues 
that the main effect of the interaction is to generate the mass of the quarks which 
become quasi-particles that interact only weakly. Neglecting this interaction 
entirely, one can view the process of mass generation as the transition of quarks 
with mass mo = to mass m resulting in a mass term mipip in the Lagrangian 
(301). We thus find that the dynamically generated mass is determined by the 
gap equation 

m = -2G(0|^V|0) • (302) 

How do we actually calculate the condensate (O|^0V|O)? To this end we go back to 
(296) and express the condensate in terms of the full propagator 'at the origin', 
i.e. at space-time point x — 0, 

(0|^V|0) m = -itvS F (x = 0) = -itr / -^j- 1 —— . (303) 

J (2TTpp — m + ie 

The following remarks are in order. First wc note that within mean-field ap- 
proximation the dynamically generated mass m is constant, i.e. independent of 
the momentum p. Furthermore, the full propagator Sf is defined in terms of 
the constituent mass m, so that the gap equation (302) becomes an implicit 
self-consistency condition where the mass m to be determined appears on both 
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sides. This equation will be solved in a moment. The integral appearing on the 
r.h.s. of (303) is quadratically divergent 14 and has to be regulated. If we use a 
cutoff A we have, for dimensional reasons, 

(0|V^|0) ~ A 2 m . (304) 

A particularly intuitive way to calculate the condensate (which will also be 
used in the light-cone case) is based on the Hellmann-Feynman theorem. This 
states that 'the derivative of an expectation value is the expectation value of the 
derivative', if the expectation is taken between normalizable states. If applied to 
the vacuum expectation value of our mean-field Hamiltonian, 

H{m) =Hq + rmptp , (305) 

which is the vacuum energy density £{m) = (0|W(ra)|0), the theorem yields, 

(0|^|0) ro = ±S(m) = - A / |f,2 v ^T^ . (306) 

The factor of two in the integrand is due to the fermion spin degeneracy. If we 
choose a noncovariant three-vector cutoff, |p| < A 3 , the result for the condensate 
is of the form (304), 



(<W|0>, 



mA 



2ir 2 



2 



1 + S - In S 2 + 0(5 2 ) 



(307) 



with S = m/A^. The dynamical mass is found by inserting this result into the 
gap equation (302). A nontrivial solution m = m(G) ^ arises above a critical 
coupling G c which is determined by the identity m(G c ) = 0. It is known from 
the theory of critical phenomena that the mean-field approximation leads to a 
square root behaviour of the mass around G c = An 2 / A 2 , 

m{G) ~ (G - G c ) 1/2 , G > G c . (308) 

We thus have seen that the NJL model describes the transformation of bare 
quarks of mass mo(= 0) into dressed (or constituent) quarks (Q) of mass m =/= 
0. By considering the bound state equation in the pseudoscalar channel one 
can also verify Goldstone's mechanism: there is a massless QQ bound state, 
which is identified with the pion, exactly if the gap equation holds (Nambu and 
Jona-Lasinio 1961; Vogl and Weise 1991). This concludes the presentation of 
spontaneous chiral symmetry breaking within the model. One should keep in 
mind that it follows the same pattern as in QCD. 

Let me conclude the NJL 'crash course' by reemphasizing that the model does 
not confine. This means in particular, that there is a nonvanishing probability 
for mesons to decay into constituent quarks. We thus cannot expect to obtain 



14 As the condensate involves the product of field operators at coinciding space-time 
points, this clearly is a short-distance singularity. 
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reliable estimates for strong decay widths of mesons or other quantities that are 
not dominated by their chiral properties. 

It is getting time to discuss the light-cone formulation of the model. For 
a 'light-cone physicist' the model is interesting for several reasons. I have al- 
ready pointed out that the lack of renormalizability is welcome because we only 
have to worry about a proper regularization. Furthermore, the model addresses 
the conceptually important questions of spontaneous symmetry breaking and 
condensates. Are these in conflict with a trivial vacuum? Finally, a constituent 
picture seems to be realized which should make a truncation of the light-cone 
Schrodinger equation feasible. Now, we again encounter the standard rule of 
physics that there is nothing like free lunch. Here, this is mainly due to the ap- 
pearance of complicated constraints for part of the fermionic degrees of freedom. 
Let me thus make a small aside on the special features of light-cone fermions. 

Light— Cone Fermions. The solution of the Dirac equation (for free fermions 
of mass mo) has the following light-cone Fock expansion at x + = 0, 

4>(x,0) = J2l°°j£- I [b(k,\)u(k,\)e- ik - x +tf(k,\)v(k,\)e ik - x 

(309) 

where we recall the notations 

k = (k + ,kj_) , je=(x~,xj_), k ■ x = ^k + x~ — k± ■ x± . (310) 

Like for scalars, the Fock measure is independent of the mass m . The Fock 
operators satisfy the canonical anti-commutation relations, 

{b(k, A) , tf(p, A')} = {d(k, A) , rft(p, A')} = 167r 3 fc + <5 3 (fc - p) , (311) 

The basis spinors u and v obey the Dirac equations, 

(#-m )u(fc,A) = 0, (312) 
(# + m )v(fc,A) = 0, (313) 

and are explicitly given by 

u(fc,A) = -^{k + +(3m + a i k i )X x , (314) 
Vk+ 

v{k,X) = -^{k+ - j3m + aiki)X- X ■ (315) 
vk + 



The four-spinor X will be defined in a moment; on and (3 are the standard 
hermitean Dirac matrices. The crucial point is the decomposition of the fermion 
field, ip = i[> + + tp- into 'good' (+) and 'bad' (— ) components, tp± = A±ip, by 
means of the projection matrices, 
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The spinor X appearing in (314) and (315) is an eigenspinor of A + , A + X\ = X\. 
The Dirac equation decomposes accordingly into two equations. The one for Th- 
reads 

2id+tp- = (-ifd 1 + m )7 + ^+ , (317) 

and does not contain a light-cone time derivative. Therefore, the bad component 
ip- is constrained; it can be expressed in terms of the good component, i.e. ip- = 
ip-[il) + ]. Again, this requires the inversion of the notorious spatial derivative d + . 
As a result, the free Dirac Hamiltonian (density) only depends on ?/>+, 

n = V-t ~^ + m§ ^+ = w hM > ( 318 ) 

where we easily recognize the light-cone energy, {k\ + m^)/k + . 

It turns out that in case of a four-fermion interaction like in the NJL model 
the constraint becomes rather awkward to solve. In particular, its solution has 
to be consistent with the mean field approximation employed (Dietmaier et al. 
1989). This can be achieved most elegantly by using the large- N expansion 
(Bentz ct al. 1999; Itakura and Maedan 2000). Nevertheless, the light-cone 
Hamiltonian of the model is a rather complicated expression . We will there- 
fore follow an alternative road which is the topic of the next subsection. 

5.3 Schwinger Dyson Approach 

The first derivation, analysis and solution of a light-cone bound-state equation 
appeared in 't Hooft's original paper on what is now called the 't Hooft model 
('t Hooft 1974). We have discussed this model in the last section where we also 
rederived 't Hooft's solution. Interestingly, 't Hooft did not use the light-cone 
formalism in the manner we presented it and which nowadays might be called 
standard. This amounts to deriving the canonical light- cone Hamiltonian and 
setting up the associated system of bound-state equations by projecting on the 
different sectors of Fock space (cf. Section 4). Instead, he started from covariant 
equations, namely the Schwinger-Dyson equations for the quark propagator (or 
self-energy), and the Bethe-Salpeter equation for the bound-state amplitude, 
which needs the quark self-energy as an input. The light-cone Schrodinger equa- 
tion was then obtained by projecting the Bethe-Salpeter equation onto hyper- 
surfaces of equal light-cone time. In this way, one avoids to explicitly derive the 
light-cone Hamiltonian, which, as explained above, can be a tedious enterprise 
in view of complicated constraints one has to solve. Let us therefore have a closer 
look at this way of proceeding. 

The Schwinger Dyson Equation for the Propagator. The first step in 
the program 15 is to solve the Schwinger-Dyson equation for the propagator, 

15 For recent literature on the Schwinger-Dyson approach, see e.g. (Roberts and 
Williams 1994; Roberts and Schmidt 2000; Maris and Tandy 2000). 
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or, cquivalcntly, for the quark self-energy. As this cannot be done exactly, one 
resorts to mean-field (or large-TV) approximation. This is essentially what has 
been done in the last subsection. Let us rewrite this in terms of Schwinger-Dyson 
equations. The one for the full propagator S reads 

S = S + S SS , (319) 

and is formally solved by 

S(P) = q - 1( * w : , (320) 
where So is the free propagator, 

p-m 

and S the quark self-energy. In mean-field approximation, it is momentum 
independent and defines the constituent mass through the gap equation, E = 
const = m, see (302) and (303). 

To solve the latter in the light-cone framework, we basically just have to 
calculate the condensate. As in the standard approach, this can be obtained 
via the Feynman-Hellman theorem by differentiating the energy density of the 
quasi-particle Dirac sea, 

— oo 

oo 

o 

Again, as it stands, the integral is divergent and requires regularization. In the 
most straightforward manner one might choose m 2 / A < k + < A and |kj_| < A, 
so that the condensate becomes 

A A 2 

<0|lM0> m = -£ / t?/^i) = -^ 21 ^- (323) 

m 2 /A 

Plugging this result into the gap equation (302) one finds for the dynamical mass 
squared, 

m 2 (G)=^ 2 exp(^-U) . (324) 

The critical coupling is determined by the vanishing of this mass, m(G c ) = 0, 
and from (324) we find the surprising result 

G c = . (325) 
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This result, however, is wrong since one knows from the conventional treatment 
of the model that the critical coupling is finite of the order 7r 2 /A 2 , both for 
covariant and noncovariant cutoff (Nambu and Jona-Lasinio 1961). In addition, 
it is quite generally clear that in the free theory (G = 0) chiral symmetry is 
not broken (as m = 0) and, therefore, this should not happen for arbitrarily 
small coupling, either, cf. (308). The remedy is to use the invariant-mass cutoff 
(Lepage and Brodsky 1980), 

M » s !^y^' < 326 » 

where we have denned the longitudinal momentum fraction, x = k + / A. This 
provides a cutoff both in x (or k + ) and k±, 

< fcl < A 2 x{l - x) - m 2 , (327) 
x < x < xi , (328) 
x ,i = |(1 T Vl-4e 2 ) , (329) 

with e 2 = m 2 / A 2 . Note that the transverse cutoff becomes a polynomial in x. 
The fc_L-integration thus has to be performed before the x-integration. 

For the condensate (322) the invariant-mass cutoff (326) yields an analytic 
structure different from (323), 



777/12 

(0|^0} ro = --^r(l + 2e 2 me 2 + O(e 2 )) , (330) 



where we have neglected sub- leading terms in e 2 . The result (330) coincides with 
the standard one, (307), if one identifies the noncovariant cutoffs according to 

A 2 = 4(A 2 3 + m 2 ) . (331) 

This has independently been observed by (Bentz et al. 1999). From (330), one 
infers the correct cutoff dependence of the critical coupling, 

4tt 2 

G c = ■ (332) 

The moral of this calculation is that even in a nonrenormalizable theory like the 
NJL model, the light-cone regularization prescription is a subtle issue. 

In the NJL model with its second-order phase transition of mean-field type, 
the usual analogy with magnetic systems can be made. Chiral symmetry cor- 
responds to rotational symmetry, the vacuum energy density to the Gibbs free 
energy, and the mass m to an external magnetic field. The order parameter 
measuring the rotational symmetry breaking is the magnetization. It is obtained 
by differentiating the free energy with respect to the external field. This is the 
analogue of expression (322) derived from the Feynman-Hellmann theorem. 
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The Bound— State Equation. Once the physical fermion mass m is known by 
solving the gap equation, it can be plugged into the Bethe-Salpeter equation for 
quark-antiquark bound states (mesons), given by 

Xbs = SiS 2 Kxbs • (333) 

Si, S2 denote the full propagators of quark and anti-quark, K the Bethe- 
Salpeter kernel, and xbs the Bethe-Salpeter amplitude. From the latter, one 
obtains the light-cone wave function via integration over the energy variable k~ 
(Michael and Payne 1982; Brodsky et al. 1985; Chakrabarty et al. 1989; Liu and 
Soper 1993), 

— X Bs(fc), fc=(fc + ,k ± ). (334) 

In ladder approximation (again equivalent to the large-iV limit), (333) and 
(334) become 



{k) = J ^ S ^ k P ) / (03 / %K{k,p) X Bs(p) , (335) 



with P denoting the bound-state four- momentum. On the left -hand-side, the 
projection onto x + — (i.e. the fc~-integration) has already been carried out. 
On the right-hand-side, the two integrations over k" and p~ still have to be 
performed. Whether this can easily be done depends of course crucially on the 
kernel K, which in principle is a function of both energy variables. For the NJL 
model, however, K assumes the very simple form, 

K(k,p) = 2 75 ® 75 -7^75 7^75, (336) 
i.e. it is momentum independent due to the four-point contact interaction, 

W~ J d 4 x J d i y{i>r^){x)5 i {x-y){i)ri;){y) . (337) 

Thus, the p~ -integration immediately yields (f>LC, and the k~ -integration can be 
performed via residue techniques and is completely determined by the poles of 
the propagators, S(k) and S(k — P). As a result, one finds a nonvanishing result 
only if < k + < P + , and one of the two particles is put on-shell, e.g. k 2 = m 2 , 
as already observed by (Gross 1989). 

The upshot of all this is nothing but the light-cone bound-state equation, 
which explicitly reads 

^ LG(x ' k±) = VP-M 2 

x J d V J -£^ tr [750Lc(j/,P_L)] A (y,p ± ) 

G $ + mh^ 5 $-f+m) 
x(l -x) M 2 — M 2 



Light-Cone Quantization: Foundations and Applications 71 

(338) 

Here we have defined the longitudinal momentum fractions x = k + /P + , y = 
p + /P + , the on-shell momentum k = (fc~,kj_, k + ) with k~ = {k\_ + m 2 )/k + , and 
the bound state mass squared, M 2 = P 2 , which is the eigenvalue to be solved 
for. 6a{x,Pj_) denotes the invariant mass cutoff (326). 

Now, while (338) may appear somewhat complicated it is actually very sim- 
ple; indeed, it is basically already the solution of the problem. The crucial ob- 
servation is to note that the two integral expressions are mere normalization 
constants, 

f d 2 

Ca = J dy J -g^ 1 tr[7 5 (/> L c(j/,Pj.)]0yi(j/,Pj.) , (339) 

D A P^j\y J ^tr[7^ L cO/,pJ]My>Pj ■ (340) 
Thus, the solution of the light-cone bound-state equation (338) is 
, , . , 2GC A (^ + m) 75 (^-f+m) 

<PLC(X,k ± ) = —2 —2 

x(l — x) M z — Mq 

GD A $ + m)flS-f+™) mi x 
+ x(l-x) M 2 -M 2 ' [ ' 

with yet undetermined normalization constants Ca and D A - As a first check 
of our bound-state wave function (34f ) we look for a massless pion in the chi- 
ral limit. To this end we decompose the light-cone wave function into Dirac 
components according to (Lucha et al. 1991), 

4>LC = <f>s + 0p75 + <7m75 + ^7m + #f <V • ( 342 ) 
Multiplying (341) with 75, taking the trace and integrating over k we find 

_ GC A f 1 , f , 2 , M 2 x + M 2 (l-x) 

Ca = -^J J x (l — x)(M 2 — Mq) X 

+ iS^ 2 J dx J <? k ± e ^ x ^) ■ ( 343 ) 

In the chiral limit one expects a solution for M = 0, the Goldstone pion. In this 
case one obtains 

1 = 2^ J dx J d2k ^ x ' • ( 344 ) 

This is exactly the gap equation (302) using the definition (322) of the conden- 
sate (with the invariant-mass cutoff (326) understood in both identities). Note 
once more the light-cone peculiarity that the (Fock) measure in (344) is entirely 
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mass independent. All the mass dependence, therefore, has to come from the 
(invariant-mass) cutoff. Otherwise one will get a wrong behavior of the dynam- 
ical mass m as a function of the coupling G, as was the case in (324). 

With this in mind, we see that the Goldstone pion is a solution of the light- 
cone bound-state equation exactly if the gap equation holds. This provides ad- 
ditional evidence for the self-consistency of the procedure. The deeper reason 
for the fact that the quark self-energy and the bound-state amplitude satisfy 
essentially the same equation, is the chiral Ward identity relating the quark 
propagator and the pseudoscalar vertex (§avkh and Tabakin 1998). 

Our next task is to actually evaluate the solution (341) of the bound-state 
equation. 0lc is a Dirac matrix and therefore is not yet a light-cone wave func- 
tion as defined in Section 4. The relation between the two quantities has been 
given by (Liu and Soper 1993), 

2P + ^(x,k ± ,A,A') =u(xP + ,k ± ,\) 1 +fac(kh + v(xP\~k ± ,\') , (345) 

where we have denoted x = 1 — x to save space. A somewhat lengthy calculation 
yields the result (Heinzl 1998), 

2QP+ /\/xx~ I ' 2C 'a \ 
4>(x, k ± , A, A') = M2 _ M 2 {-jfU X M lb v x , - D A u x f m , J , (346) 

with the arguments of the spinors u and v suppressed. At this point we have to 
invoke another symmetry principle. (Ji et al. 1992) have pointed out that the 
spin structure (w-Tv) should be consistent with the one obtained form the instant 
form spinors via a subsequent application of a Melosh transformation (Melosh 
1974) and a boost. Using this recipe, one obtains the following relation between 
the constants Ca and Da, 

2C A /M = -D A = N/2G . (347) 

As a result, the spin structure in (346) coincides with the standard one used 
e.g. in (Dziembowski 1988; Chung ct al. 1988; Jaus 1990; Ji and Cotanch 1990; 
Ji et al. 1992). The NJL wave function of the pion thus becomes 



N P + I \/tt 

VOr, kj_, A, A') = M2 ^ u x (M + f) l5 v x , 9(A 2 - M 2 ) . ; 34<S) 



Not surprisingly, the off-shcllncss M 2 — M 2 appears in the denominator. N is 
the normalization parameter defined in (347), and the spin (or hclicity) structure 
is given by 

u(xP+, k ± , A) (M + fi 75 v(xP\ -k ± , A') = 

=L — \\(mM + m 2 - k 2 ± + M 2 xx)S x _ v - k_ x (M + 2m)8 xx ] (349) 
xxP+ L V / J 

where we have used (314), (315) and denoted k x = k 1 +i\k 2 . The first term with 
spins anti-parallel corresponds to L z = 0, the second one (with spins parallel) to 
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L z = ±1. It has already been pointed out by Leutwyler that both spin alignments 
should contribute to the pion wave function (Leutwyler 1974a; Leutwyler 1974b). 
Note that the latter is a cutoff dependent quantity. This is necessary in order to 
render the wave function normalizable. A single power of the off-shellness in the 
denominator is not sufficient for that. Only the cutoff guarantees the boundary 
conditions (251) so that the wave function drops off sufficiently fast in x and k±. 

As we are interested in analyzing the quality of a constituent picture, we 
approximate the pion by its valence state, denoting ip = tp2, 

\ 7T:P )=J2j Q dx J Mx, k±, A, A') \qq : x, kj_, A, A') (350) 

which should be compared with the general expression (240) . The normalization 
of this state is given by (237), 

(tt : P'\tt : P) = 16ir 3 P+5 3 {P - P') . (351) 

As usual we work in a frame in which the total transverse momentum vanishes, 
i.e. P = (P + ,P± = 0). Expression (351) yields the normalization (247) of the 
wave function, 

/ ^hM*,k ± ,A,A')| 2 ^||?H| 2 = l. (352) 

It is of course a critical assumption that the probability to find the pion in its 
valence state is one. In this way we enforce a constituent picture by fiat, and it 
is clear that such an assumption has to be checked explicitly by comparing with 
phenomenology. 

5.4 Observables 

With the light-cone wave function at hand, we are in the position to calculate 
observables. To proceed we will employ the following two simplifications. First 
of all, we will always work in the chiral limit of vanishing quark mass, m = 0, 
which, as we have seen, leads to a massless Goldstone pion, M = 0. We write 
the pion wave function as a matrix in helicity space, 

_ * /-2m(fc'- 8 fc 2 ) m 2 -kl \ , 2 2 

(353) 

Note that, in the chiral limit, the wave function becomes independent of x (apart 
from cutoff effects) . This actually agrees with our findings in the two-dimensional 
't Hooft model, cf. (276). 
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The diagonal terms in (353) correspond to parallel spins, the off-diagonal 
ones to anti-parallel spins. The different components are related by the symme- 
try properties, 

V>2U = $2TT ' ^2iT = ~i>2M ■ (354) 

Second, we will go to the large-cutoff limit, that is, we will keep only the leading 
order in e 2 = m? / A 2 . We thus assume that the cutoff is large compared to 
the constituent mass. From the standard values, A ~ 1 GeV, m ~ 300 MeV, 
we expect that this assumption should induce an error of the order of 10%. 
The technical advantage of the large-cutoff limit is a simple analytic evaluation 
of all the integrals we will encounter. Furthermore, the leading order will be 
independent of the actual value of the constituent mass. It should be mentioned 
that the same procedure has been used in calculations based on the instanton 
model of the QCD vacuum (Petrov et al. 1999). There, the ratio e 2 can be related 
to parameters of the instanton vacuum, namely 

6 2 = {mpf ~ {p/Rf , (355) 

where p ~ 1/3 fm is the instanton size and R the mean distance between in- 
stantons. Thus, e 2 can be identified with the 'diluteness parameter' or 'packing 
fraction' of the instanton vacuum and hence is parametrically small, e 2 ~ 1/4. 

The upshot of all this is that we work with the extremely simple model wave 
function (Radyushkin 1995) 

V>2U - N 0(A 2 - M 2 ) , V 7 2TT = , (356) 

which is entirely determined by two parameters, the normalization constant N 
and the cutoff A. We thus need two constraints on the wave function to fix our 
two parameters. 



Normalization. As announced, we enforce a constituent picture by demanding 
(352) which decomposes into 

i = \\H\ 2 = ||V>2ull 2 + ll^nll 2 + IhMI 2 + Ihfcull 2 • (357) 

Explicitly, one finds 



ll*2 U || 2 = ll*21tl| 2 = « ! 



N 2 f 1 fA 2 x(l-x) N 2 A 2 

' dx / dk ± = -T^nr = 1/2 , (358) 



16tt 2 ./o Jo 96tt 2 



while the components with parallel spins have vanishing norm in the large-cutoff 
limit, ||V> 2TT || 2 = ||V>2ul| 2 =0,cr. (356). 
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Pion Decay Constant. A second constraint on the wave function is provided 
by the pion decay constant / w , which appears in the semi-leptonic process 7r — > 
\w. The relevant matrix element is 

(0|^-(0)7 + 75^(0)|7r + (P + )> = zV2P + /. , (359) 

ipa and ip u denote the field operators of the d and u quark in the pion. If we 
insert all quantum numbers, the pion state to the right of the matrix element is 
given by 



-^(KtM) _ Ki«cT)) • (360) 



The spatial (or internal) structure of the state is encoded in the light-cone wave 
function xp2 = V'Ju- If we insert the Fock expansions (309) for ipj and ip u as 
well as the pion state (350), we obtain the following constraint on the pion wave 
function, 

[o-S^^^'Wr <361) 

The left-hand-side is basically the (position space) 'wave function at the origin'. 
Quark (u) and antiquark (d) thus have to sit on top of each other in order to 
have sizable probability for decay. Note that only the L z — component con- 
tributes. Concerning the effect of higher Fock states, it can be shown (Lepage 
et al. 1981; Brodsky and Lepage 1989) that indeed only the valence wave func- 
tion contributes to (361). This constraint is therefore exact and holds beyond 
a constituent picture. Empirically, the pion decay constant is f v = 92.4 MeV 
(Holstein 1995). 

As already stated, this is our second source of phenomcnological informa- 
tion to fix cutoff and normalization. Using the explicit form (356) of the wave 
function, the constraint (361) becomes 



l! dx I 



^ , ( ^ na2 ! ^ 



With (362) and (358) we now have two equations for our two parameters which 
accordingly are determined as 

N = V3/U , (363) 
A = AttU ^ 1-16 GeV . (364) 

The value (364) for the cutoff A is the standard scale below which chiral effective 
Lagrangians are believed to make sense (Manohar and Georgi 1984). It is reas- 
suring that within our approximations we get exactly this value. This means that 
we are not doing something entirely stupid. A more severe test of consistency is 
provided by the next constraint to be satisfied. 
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Constraint from tv° — > 2j. This constraint has also been derived by Brodsky 
and Lepage (Lepage et al. 1981) within an analysis of the tt-j transition form 
factor. It assumes the very simple form, 

f dxip2n( x ^°-i-) = "T^ ■ ( 365 ) 

JO Jtt 

Inserting the light -cone wave function (356), the right-hand-side simply be- 
comes the normalization N which is indeed consistent with our findings (363) 
and (364). We mention in passing that the constraint (365) usually is the sim- 
plest way to fix the normalization N. Its derivation, however, is more complicated 
than that of (361). 



Pion Form Factor. We proceed by calculating the pion electromagnetic form- 
factor. It is defined by the matrix element of the electromagnetic current J^ m 
between pion states, 

(tt : P\J£ m \n : P') = 2(P + P'Y F(Q 2 ) , Q 2 = -(P - P') 2 . (366) 

Considering \i = + in a frame where P = (P + , 0) and P' = (P + , q_|_) one is led 
to the the Drell-Yan formula (Drell and Yan 1970; Brodsky and Lepage 1989), 

*"(<&) = £ I dx f ^V*(*,k' ± ,A,A')#r,k ± ,A,A'). (367) 

The transverse momentum of the struck quark is k' ± = k^ + (1 — x)q ± . The 
formula (367) with its overlap of two wave functions on the right-hand-side is 
rather similar to the nonrclativistic result as will be shown in what follows. 

The form factor of a nonrelativistic system is given by the Fourier transform 
of the charge distribution (normalized to one), that is, 

F(p) = J d\r{v) V>« e 4p ' x - J -0^ ?/>*(k + p) V(k) . (368) 

It is important to note that k and k' = k + p arc relative momenta, 

k= ^(^2ki -mik 2 ) , (369) 

(and analogously for k'), so that the ki.2 are the actual particle momenta. Ac- 
cordingly, p is the relative momentum transfer, 

P = k -k=— (ki-kx)^— q = a; 2 q, (370) 

where, in the last step, we have used (268). Plugging this into the form factor 
(368) we obtain the formula, 

/d 3 k 
^V>*(k + x 2 qMk), (371) 



Light-Cone Quantization: Foundations and Applications 



77 



which, as promised, is quite similar to (367). 

If one sets the momentum transfer q ± in (367) equal to zero, one is left 
with the normalization integral (352), so that the form factor is automatically 
normalized to one (the same is true in the nonrelativistic case). 

We will use the Drell-Yan formula for the form factor to determine the pion 
charge radius r v , which is given by the slope of the form factor at vanishing 
momentum transfer, 

*.2 



F(qi) = l-fqi + 0(<f ± ). 
Using (367) this results in the nice explicit formula, 



(372) 



r 2 = -- 



dx 



d 2 k ± d 2 



167T 3 dqidq 



■ipZ(k ± + xq ± ) 



V>(k ± ) 



(373) 



q±=o 



Upon inserting the wave function (356), however, one encounters a problem. 
The sharp cutoff (corresponding to a step function) is too singular to lead to a 
reasonable result. The derivatives in (373) are concentrated at the boundary of 
the support of the wave function which in the end leads to artificial infinities. 
Thus, for the time being, we resort to a smooth cutoff, 



9 s A {x,k ± ) = cxp 



kl 



A 2 x(l - x) 



(374) 



which basically transforms the sharp-cutoff wave function (356) to the Gaussian 
(286) (with i7i set to zero). Plugging this into (373) yields the pion charge radius 



r 2 = 

' IT 



12 



A 2 4tt 2 /2 



= (0.60 fm) 2 



(375) 



This is the standard result for the NJL model (Blin ct al. 1988) and has also 
been obtained within the instanton model (Diakonov 1996). It is slightly smaller 
than the experimental value, r T = 0.66 (Amendolia et al. 1986), a discrepancy 
which is usually attributed to the use of the large-cutoff limit. A pole fit using 
our value of the pion charge radius is displayed in Fig. 10. 



Transverse Size. As in Section 4, Example 3, we can use the light-cone wave 
function (356) to calculate the r.m.s. transverse momentum which leads to 

(ki) = j\x J ^k\ |hA 2U (x,k ± )|| 2 = ^(370McV) 2 . (376) 

This actually coincides with the result (288) for the smooth-cutoff Gaussian wave 
function (286) or (374). Therefore, unlike the charge radius rv, the r.m.s. trans- 
verse momentum is insensitive to the details of the cutoff procedure. We thus 
have (fc^) 1 / 2 ~ m > M w , which confirms that the pion is highly relativistic. 
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Fig. 10. The pion form factor squared vs. momentum transfer q 2 = q\. The full line 
is the monopole fit of (Amendolia et al. 1986), \F\ 2 = n/(l + q\r% /6) 2 with n = 0.991, 
r 2 = 0.431 fm 2 ; the dashed line is the same fit with our values, n = 1, r 2 = 0.36 for 2 . 
The agreement is consistent with the expected accuracy of 10%. 
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The r.m.s. transverse momentum can easily be translated into a transverse 
size scale R±, 

R\ = l/(ki) ~ (0.54 fm) 2 . (377) 

This is slightly smaller than the charge radius which we attribute to the fact that 
the charge distribution measured by the charge radius does not coincide with the 
distribution of baryon density. The 'core radius' R± is sometimes related to the 
decay constant / w via the dimensionless quantity C — f n R± (Weise 1984). In 
constituent quark models one typically gets C ~ 0.4. This implies the fairly large 
value R± ~ 0.8 fm. Using standard many-body techniques, Bernard et al. have 
calculated this quantity in a model treating the pion as a collective excitation 
of the QCD vacuum, and find C ~ 0.2 (Bernard et al. 1985). This result is close 
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to what we get from (377), 

C = f n R± ~ 0.25 . (378) 

Thus, though we work within a constituent picture, we do get a reasonable 
value. We believe that this is due to the intrinsic consistency of the light-cone 
framework with the requirements of relativity, a feature that is lacking in the 
constituent quark model. 



(Valence) Structure Function. The pion structure function arises in the 
description of deep inelastic scattering off charged pions. In terms of light-cone 
wave functions it is defined as the momentum fraction x times the sum of 'quark 
distributions' f q weighted by the quark charges e q . The quark distributions are 
given by the squares of light-cone wave functions integrated over k±. For the 
valence structure function of the pion we thus have the formula (Lepage et al. 
1981; Brodsky and Lepage 1989), 

FZ(x)=x[elfZ(x)+elr d (x)] = = / ^ |^(x, k ± , A, A')| 2 , 

A A' 

(379) 

where the /" denote the valence quark distributions. For the model wave function 
(356) the structure function becomes 

F-(x) = fx 2 (l-x), (380) 

which in turn leads to the (valence) quark distribution 

r(x) = fZ{x) = fl(x) =6x(l-x). (381) 



The following consistency checks can be made. The probability to find a valence 
quark in the pion, 

t-i 

dxf v (x) = l, (382) 

lo 

is unity, as it should. For the mean value of the momentum fraction x carried 
by one of the quarks one finds 



Jo 



(x) = / dxxf v (x) = 1/2 . (383) 
Jo 

Thus, on average, the quarks share an equal amount of longitudinal momentum, 
which again, of course, is the correct result. 

If one compares with other NJL calculations of the structure function (Shigc- 
tani et al. 1993; Bentz et al. 1999) and with the empirical parton distributions 
in the literature (Gluck et al. 1999), one finds reasonable qualitative agreement. 
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Let me finally point out that it is not entirely obvious to which actual mo- 
mentum scale our results correspond. The transverse-momentum cutoff is in- 
dependent, A 2 (x) ~ A 2 x(l — x), so, if we use the average x of (383), (x) ~ 1/2, 
a natural scale seems to be 16 



Q 



1/2 

(x)(l-(x)) A ~ A/2 ~ 600 MeV 



(384) 



Pion Distribution Amplitude. The pion distribution amplitude was orig- 
inally introduced to describe hard exclusive processes involving pions (Lepage 
and Brodsky 1979a; Lepage and Brodsky 1979b; Lepage and Brodsky 1980). The 
pion formfactor at large Q 2 , for example, is given by the following convolution 
formula, 

F(Q 2 )= f dx f dycj>*(x,Q)T H (x,y;Q)cf>(y,Q) \l + 0(1/ Q)] , (385) 



Jo Jo 

where Q denotes the large momentum transfer, Th a 'hard scattering amplitude' 
and <f> the pion distribution amplitude. While the amplitude Th is the sum of 
all perturbative contributions to the scattering process, <fi is nonperturbative in 
nature. The convolution formula is thus a prominent example where we see the 
principle of factorization into 'soft' and 'hard' physics at work. 

The pion distribution amplitude is rather straightforwardly related to the 
light-cone wave function of the pion, 

<P(x,Q) ~ J ^ e " p+z ~/>|<K0) 7 +75^-,(U)KP+)> 

^4f(^)- (386) 



/ 



The normalization is fixed by demanding that <j> integrates to unity. 

Brodsky and Lepage have shown that <j> obeys an evolution equation of the 
form (Lepage and Brodsky 1979a; Lepage and Brodsky 1980), 

Q ^Q ^ Q) = I dy V{X ' V] Q) Q) • (387) 

where the evolution kernel V is determined by perturbative QCD. For Q — > oo, 
(387) has the asymptotic solution 

</> BB (x)=6x(l-x), (388) 

(which is normalized to 1). The (nonasymptotic) pion distribution amplitude 
has been a rather controversial object. For a while people have tended to be- 
lieve in a 'double-humped' shape of the amplitude (due to a factor (1 — 2x) 2 ), 
which was originally suggested by Chcrnyak and Zhitnitsky using QCD sum rules 



I thank W. Schweiger for discussions on this point. 
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(Chernyak and Zhitnitsky 1984). In 1995, however, the CLEO collaboration has 
published data (Savinov 1995) that seemed to support an amplitude that is not 
too different from the asymptotic one (Kroll and Raulfs 1996). Theoretical ev- 
idence for this fact has recently been reported in (Belyaev and Johnson 1998a; 
Belyaev and Johnson 1998b; Petrov and Pobylitsa 1997; Petrov ct al. 1999). 
Belyaev and Johnson, for instance, have found two constraints which should be 
satisfied by the distribution amplitude (Belyaev and Johnson 1998a), 

4>{x = 0.3) = 1±0.2 , 

<t>(x = 0.5) = 1.25 ±0.25 , (389) 

which are consistent with an amplitude being close-to-asymptotic. 

Last year, the experimental developments have culminated in a direct mea- 
surement of the distribution amplitude at Fermilab (Ashery 1999). At a fairly low 
(i.e. nonasymptotic) momentum scale of Q 2 ~ 10 GeV 2 , one finds a distribution 
amplitude that is very close to the asymptotic one. 



Fig. 11. The (asymptotic) pion distribution amplitude. The vertical lines mark the 
constraints (389) of Belyaev and Johnson. 
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Let us see what we get in the NJL model. The distribution amplitude is given 

by 

<Pn.jl(x) = J V2Tl(£,k_L) = 6x(l -x) = (j) as (x) . (390) 

Thus, in the large cutoff limit, the NJL distribution amplitude is exactly given 
by the asymptotic one! Upon comparing with (381), we see that </>njl coincides 
with the quark distribution This is accidental and stems from the fact that, 
due to the use of a step function in (356), ^2 ~ IV^I 2 - 

In Figure 11 we have displayed our distribution amplitude together with the 
constraints (389), represented by the vertical lines. 

Our findings are thus consistent with the recent Fermilab experiment (Ashery 
1999). One should, however, be aware of the fact that our energy scale of Q ~ 0.6 
GeV from (384) is below the experimental one of Q ~ 3 GcV. It is also somewhat 
lower than Q ~ 1 GeV, which has been assumed by Belyaev and Johnson in 
their analysis of the distribution amplitude in terms of light-cone quark models 
(Belyaev and Johnson 1998b). 

6 Conclusions 

In this lecture I have discussed an alternative approach to relativistic (quantum) 
physics based on Dirac's front form of dynamics. It makes use of the fact that for 
relativistic systems the choice of the time parameter is not unique. Our particular 
choice uses null-planes tangent to the light-cone as hypersurfaces of equal time. 
This apparently trivial change of coordinates has far-reaching consequences: 

— The number of kincmatical (i.e. interaction independent) Poincare genera- 
tors becomes maximal; there are seven of them instead of the usual six, 
among them the boosts. 

— Lorcntz boosts in z-direction become diagonal; the light-cone time and 
space coordinates, x + and x~ , respectively, do not get mixed but rather 
get rescaled. 

— As a consequence, for many-particle systems one can introduce frame-indep- 
endent relative coordinates, the longitudinal momentum fractions, Xi, and 
the relative transverse momenta, kj_j. 

— Because of a two-dimensional Galilei invariance, relative and center-of-mass 
motion separate. As a result, many formulae are reminiscent of nonrelativis- 
tic physics and thus very intuitive. 

— This is particularly true for light-cone wave functions which, due to the last 
two properties, arc boost invariant and do not depend on the total momen- 
tum of the bound state. They are therefore ideal tools to study relativistic 
particle systems. 

— The last statement even holds for relativistic quantum field theory where one 
combines the unique properties of light-cone quantization with a Fock space 
picture. The central feature making this a reasonable idea is the triviality 
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of the light-cone vacuum which accordingly is an eigenstate of the fully 
interacting Hamiltonian. This implies that the Fock operators create the 
physical particles from the ground state. 

As expected, however, the principle of conservation of difficulties is at work 
so that there are problems to overcome. Of particular concern to us was the 
'vacuum problem'. In instant-form quantum field theory, especially in QCD, 
many nonperturbative phenomena are attributed to the nontriviality of the vac- 
uum which shows up via the appearance of condensates. These suggest that 
the instant-form vacuum is a complicated many-body state (like, e.g. the BCS 
ground state). In addition, many of these condensates signal the spontaneous (or 
anomalous) breakdown of a symmetry. The conceptual problem which arises at 
this point is to reconcile the existence of condensates with the triviality of the 
light-cone vacuum. 

The idea put forward in these lectures is to reconstruct ground state proper- 
ties from the particle spectrum. The latter is obtained by solving the light-cone 
Schrodinger equation for masses and wave functions of the associated bound 
states. For a relativistic quantum field theory, this, in principle, amounts to 
solving an infinite system of coupled integral equations for the amplitudes to 
find an ever increasing number of constituents in the bound state. Experience, 
however, shows that the light-cone amplitudes to find more than the valence 
quanta in the bound state tend to become rather small. Note that the same is 
not true within ordinary, that is, instant-form quantization. 

A first explorative step towards solving a realistic light-cone Schrodinger 
equation was performed using an effective field theory, the NJL model. We have 
seen that, though we made a number of approximations, in particular by en- 
forcing a constituent picture, a number of pionic observables are predicted with 
reasonable accuracy, among them the pion electromagnetic form factor, the pion 
charge and core radius and the pion valence structure function at low normaliza- 
tion scale. The pion distribution amplitude (in the chiral and large-cutoff limit) 
turns out to be asymptotic. 

The results presented in these lectures provide some confidence that, also 
for real QCD, light-cone quantization may provide a road towards a reasonable 
constituent picture, in which hadrons are consistently described as bound states 
of a minimal number of constituents. How this hope can be turned into fact 
remains to be seen. 
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